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Abstract. We construct the upper and lower Bellman functions for the U" (quasi)-norms 
of BMO functions. These appear as solutions to a series of Monge- Ampere boundary value 
problems on a non-convex plane domain. The knowledge of the Bellman functions leads to 
sharp constants in inequalities relating average oscillations of BMO functions and various 
BMO norms. 



1. Introduction 

For a measurable set E C with finite, non-zero Lebesgue measure \E\, and a locally 
integrable real-valued function ip, let (9?)^ denote the average of (p over E, 

For a function 93, a cube Q, and p > 0, let ||v3||lp(q) denote the Q-normalized (quasi-)norm 
of ^p, 

Observe that ||<^||lp(q) is increasing in p : 

(1-1) IIv'IIlpi(q) < II¥'IIlp2(q) 

for pi < P2, with equality happening if and only ii ip = const on Q. 
Fix a cube Q in and let BMOP(Q) be the (factor-)space 

(1.2) BMOP{Q) = {ip£L^{Q): {\ip - {^) j\p) j < < 00, V cube J C Q}, 

with the smallest such C being the corresponding norm (quasi- norm for < p < 1 ) , 

(1-3) IIv'IIbmopcq) = sup {\(p - {(f) . 

cube JcQ 



It is known that all p-based norms defined by (1.3) are equivalent and so (1.2) defines the 
same space for all p > 0. This fact is usually seen as a consequence of the John-Nirenberg 
inequality, although using that inequality to prove it will produce suboptimal constants of 
norm equivalence. One of the primary motivations of this work is to quantify this equivalence 
precisely, in dimension 1. To this end, we relate all BMO^ norms to the BMO^ norm. The 
reason BMO^ norm plays a central role here is that it allows us to take advantage of the 
self-duality of L?{Q). From now on, we reserve the name BMO for BMO^ : 

BMO(Q) = {(^ G L\Q) : {ip^) j - (99)2 < < 00, V cube J C Q}. 

Thus, we would like to find the best constants Cp, Cp in the double inequalities 

(1.4) Cp 1 1 (/J 1 1 BMO(Q) < sup {W - {'^) j\'^)^J^ < Cp\\<p\\jiMO{Q) , P>0. 

interval JcQ 
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Some cases are trivial: (1.1) implies that Cp = 1 for p < 2, and Cp = 1 for p > 2 (the 
equalities hold, for instance, for any if with zero average and constant modulus on Q). 

To find Cp and Cp, we estimate, for any ip G BMO(Q), the quantity {\ip — {^)q\^)q (we 

will refer to this as the p-oscillation of f over Q) in terms of its 2-oscillation (v'^)r) ~ (v)^ 

and ||95||bmo{Q)- Out estimates are sharp for any p > and one can use them to estimate the 
p-oscillation of a function in terms of its g-oscillation and BMO norm, for any < p,q < oo. 
What is surprising is that, due to the nature of the optimizers in the p o 2 inequalities, we 
obtain sharp p q inequalities, whenever p E [1,2] and q S [2,co). 



The oscillation estimates immediately yield the norm equivalence statement (1.4). However, 
the norm equivalence itself may or may not be sharp; at this point we can only show sharpness 
for p > 2. 

The principal step in getting oscillation estimates is to obtain sharp two-sided inequalities 
for the p-th. power of L^(Q)-norms of BMO functions, i.e. expressions of the form 
It turns out that these quantities are always finite, meaning that BMO(Q) C LP{Q), as sets. 
These sharp estimates on BMO are our main goal. 

Following the template of the John-Nirenberg project [SVJ, we define the upper and lower 
Bellman functions for the problem: for p > and e > 0, let 



(1.5) B,,pix)= sup {{\ipn : {^) =xi,{^^)^ 



X2 



(1.6) 6,,p(x) = inf I (Iv'Dq : {v>)q = ^i' (¥'')q = ^2} , 

where BMOe(Q) is the e-baU in BMO(Q), 

BMO,(Q) = {ipG BMO(Q): ||<^||bmo(q) < e}. 

It easy to check that these functions are independent of the interval Q. Their domain is 

= {x : xf < X2 < Xi + e^}. 

Indeed, for every ip £ JiMO^{Q) and every subinterval J of Q the corresponding Bellman 
point {{if) J , {if^) j) is in : the first inequality is just Holder's inequality and the second 
one holds since ||<^||B]yjo(Q) < £• This is the same domain as in jSVj . and, as was the case 
there, each Bellman function will satisfy a Monge-Ampere equation on 0^. However, unlike 



the Bellman functions in [SV], the functions defined by (1.5) and (1.6) do not have additive 
homogeneity, while the domain fig does not allow for multiplicative homogeneity. Thus we 
cannot hope to reduce the partial differential equations to ordinary ones, as was done in the 
previous work. Instead, we construct special foliations of il^ by straight-line characteristics 
along which each Bellman function must be linear. These foliations, different for various ranges 
of p, allow us to solve the PDE, thus obtaining the Bellman functions, and, simultaneously, 
construct optimizers in the inequalities being proved. 

Having introduced the main objects of study, let us say a bit about the method. The term 
"Bellman functions" alludes to similar extremal constructs in the dynamic programming of 
R. Bellman |Bej . There are deep parallels between the Bellman frameworks of optimal sto- 
chastic control and harmonic analysis (see, for example, |NTV3j and |Vol] ) . although we make 
no use of those connections here. In the mid- to late 1980s, R. Burkholder ( |BH lB2]) started 
using specially designed functions with delicate size and concavity properties to prove sharp 
inequalities for martingales by induction on scales. In the 1990s, the method took its modern 
form in the work of F. Nazarov, S. Treil, and A. Volberg, starting with [NTVlj . |NT| . and 
|NTV2| . However, the first exact Bellman functions, as explicit solutions of extremal prob- 



lems such as (1.5) and (1.6), did not appear until jVlj. Bellman analysis on BMO originated 



with |SV| and continued in [y2] . Other notable explicit Bellman functions appeared in the work 
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of A. Melas and co-authors (see, for instance, [Mil IM2l IMNj ). although those functions were 
found using combinatorial analysis of the operator in question, the dyadic maximal function, 
as opposed to solving the Bellman partial differential equation. The Bellman functions for the 
maximal operator were taken up again in [SSt ] and [SSVj. Those works, together with |VV2j . 
initiated the study of the relationship among explicit Bellman functions, the Monge-Ampere 
geometry of Bellman domains, and the structure of optimizers in the corresponding inequal- 
ities. Most recently, that line of investigation continued in |VV1| . The current work also 
presents that unified vision. 



How do we proceed from estimates on {\^\^)q to the norm equivalence statement (1.4)? By 
the definitions of and ft^, we have, for every ip £ BM0e((5) ^-^id every subinterval / of Q, 

Replacing (p with (p — {'p) ^ gives sharp oscillation estimates: 

Analyzing each inequality separately, we get the desired norm inequalities: 

hv\lp{^M\)< sup < S|j<^|j_p(0, 

interval /CQ 

where || • || =^ || • ||bmo(Q)- 

The functions B^^p and h^^p can be viewed as special cases in a more general framework. 
Namely, take a function / on M and define, formally, the Bellman functions 

(1.7) B^j{x)= sup {{f{Lp))^: {lp)^=xi,{'p^)j=X2] 
and 

(1.8) bej{x) = inf {(/((^)), : {^), = XI, = X2} . 

(Here we explicitly allow for the possibility that one or both of these functions take on infinite 
values.) We will see that such a general view is beneficial: we will develop several canonical 
building blocks, each defined in a sub-domain of these blocks, when appropriately arranged 



and glued, produce the functions (1.7), (1.8) for various choices of /, including the power 
function. Equally important, the optimizers in the Bellman functions — those functions ip on 
which supremum or infimum is attained — turn out to be determined locally by the canonical 
blocks, rather than by what specific / is being considered. In a very tangible sense, many 
BMO inequalities have the same optimizers. 



Let us outline several important choices of / : f{s) = |s|^ yields definitions (1.5) and (1.6) 



/(s) = produces the Bellman setup for the integral John-Nirenberg inequality from [S V] : 



since limp-^o{\ip\P)^P = e ^'"^'^Q^ the choice /(s) = log |s| gives the appropriate limiting setup 



for (1.5) and (1.6), although we will see that the resulting "BMO*^," is not, in fact, BMO. 

There are at least two more choices of importance: /(s) = X(-oo,-A]u[a,oo)('5) for A > and 
/(s) = e'^l. The former yields the setup for the classical, weak-form John-Nirenberg inequality 
and the latter, for the two-sided integral John-Nirenberg inequality. Each of these two cases 
requires a slight modification of our building blocks. The first was considered in [V2j, in a 
much more limited context, and the second will be considered elsewhere. 

What are the assumed conditions on /? Since our focus in the present work is on two-sided 
inequalities, we will assume that / is even. In order to simplify exposition we will also assume 
that / is non-negative and smooth, except, possibly, at 0. 

Although all definitions in this section are valid in any dimension, at present we are only 
able to find the Bellman functions in the one-dimensional case, where all cubes in the definition 
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of BMO are intervals. We can obtain meaningful dimension-dependent estimates of the norm- 
equivalence constants in higher dimensions, but it is apart from our main interest here, which 
is in sharpness and explicit Bellman functions. In fact, we hold out the possibility that the 
Bellman functions — and so the norm estimates — are dimension-free, which would have major 
implications for analysis on BMO and related function classes (such as Ap). At this time, we 
cannot show it, since a key geometric ingredient in our proofs works only in dimension 1; we 
hope to be able to give a definitive answer for higher dimensions in the future. 

The geometry of Monge-Ampere foliations plays a central role in the construction of Bellman 
functions and their optimizers. As such, it is given a central role in our exposition. The picture 
of a foliation concisely captures the nature of the extremal problem at hand, and it is those 
foliations, first built locally and then carefully glued together, that we would like the reader 
to remember. While we provide precise algebraic descriptions of the Bellman functions, our 
proofs often appeal to their geometric nature. For example, when building optimizers, it is 
certainly possible to show that a given function is in BMO^ by a direct calculation. However, 
it is more geometrically meaningful — and often simpler — to show that all of its Bellman 
points are in fi^. 

Lastly, proving sharp norm estimates, while a lofty goal, does not require one to know the 
origin of the Bellman functions or of their optimizers; a relatively straightforward verification 
that these are, in fact, optimal would suffice. However, we choose to present their construction 
in considerable detail. This serves our second major goal: to provide a complete account of the 
modern Bellman-Monge-Ampere approach to problems with non-convex Bellman domains. 
Many other traditional harmonic analysis questions give rise to such domains and much of 
what follows should be applicable there. 

The paper is organized as follows: in Section[2j we state the sharp inequalities proving which 
was our main motivation; in Section [3j we derive the boundary value problems our Bellman 
functions should solve and outline the geometric approach to finding the solutions; in Section |4] 
we construct local Bellman candidates for several subdomains of fig; these are glued together 
to produce global solutions in Section [5j in Sections [6] and [7] we use induction on scales and 
optimizers to prove that the global candidates are, in fact, the true Bellman functions; in 
Section [8] we provide the proofs of the inequalities from Section [2| finally, in Section [9j we 



briefly consider several additional choices of / in (1.7) and (|l.8|) and state the corresponding 
Bellman functions. 

2. Sharp inequalities 



Although the main results of this work are the explicit expressions for the functions (1.5) 
and (1.6) for various ranges of p, these are relatively complicated and stated in Section[5| after 



the relevant notation has been introduced and various blocks that comprise these expressions 
have been developed. In this section, we state the immediate consequences of our knowing the 
Bellman functions: the appropriate sharp inequalities for BMO. Since we know both the upper 
and lower functions, we naturally get two-sided inequalities; even though some are elementary, 
they are included for symmetry and to emphasize their source. 

Theorem 2.1. For an interval Q and any ip € BM0((5) such that \\^p\\bmo(Q) 7^ we have: 
^/0 < p < 1, 2P-^^r^-^^^^^i{^^)^ - (^)p < (Iv. - (c,)^^)^ < - {^)lf"; 

^/l < p < 2, ir(p)||c,|r3-Jo(Q)((,.2^Q _ ^^^p < _ < (^^2^^ _ ^^^2j/\^ 

^f2<p<oo, - < (1^ - < i mMluoiQ)ii^')Q - (V')g)' 
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and these inequalities are sharp and attainable. 

As a straightforward corollary of this theorem we obtain inequalities relating BMO^ and 
BMO^ norms; for p > 2 these are the best possible. 

Theorem 2.2. For an interval Q and if G BMO(Q), we have: 

ifO<p<l, 2^-^/P\\ip\\-BMO{Q) < II'/5||bmop(Q) < II'/^IIbmo(Q); 

ifl<P<2, (f r(p))^''^ ||v?||bmo(Q) < llv'llBMOfCQ) < II'/'IIbmocq); 

if2<p<oo, \\(p\\bmo{Q) < IIv5||bmop(Q) < (i r(p))^'^^ ||(^||bmo(Q)- 
The right-hand side inequalities for p < 2 and both left- and right-hand side inequalities for 
p > 2 are sharp and attainable. 

Remark 2.3. At this point we do not know if the left-hand inequalities for p < 2 are sharp. 



Of course, Theorem 2.2 allows us to relate different BMO^ norms bypassing BMO^, but 



the resulting inequalities may no longer be sharp. However, if we instead use Theorem 2.1 and 
our knowledge of the optimizers in its statements, we get the best possible inequalities relating 
p-oscillations for certain ranges of the parameter p. Specifically, we have 



Theorem 2.4. Fix an interval Q and numbers pi € [1,2] and p2 G [2,oo). Then, for any 
(f £ BMO(Q), we have 

PiT{pi 

and these inequalities are sharp and attainable. 



As mentioned in the introduction, in this paper we do not attempt to prove sharp two-sided 
John-Nirenberg inequality, that is we do not find the Bellman functions (1.7) and (1.8) for 

f{s) = e'^l. However, because e''^' = 1 + |s| -|- X]fcl2 ^ upper Bellman functions Bp^^ 

for p > 2 turn out to have the same optimizers, we can sharply estimate 

Coupling the resulting estimate with the inequality {\ip — (ip) q\) ^ < ||v'||bmo(Q)5 'we obtain 
the following "almost sharp," two-sided integral John-Nirenberg inequality: 

Theorem 2.5. Take an interval Q and let f € BM0((5) be such that \\ip\\ = \\^\\bmo(Q) < 1- 
Then we have 

(el^-(^)Ql)^<C(||^||). 
Here the bound 1 on ||9j||bmo is sharp and the best (smallest) value of C(||99||) satisfies 

1 _ M 1 _ MP 

(2.1) < C(ll^ll) < 



1 — \\ip\\ 1 — \\ip\\ 

By "sharp" we mean that there exist functions 93 with norm 1 for which the inequality 
fails. The reader can compare this result with the one in |SV| . where the sharp one-sided 
John-Nirenberg inequality was proved. One can get a sub-optimal two-sided inequality from 
a one-sided one by simply doubling the constant. If we do that with the result in [SV| . we get 

cm)< 



which is worse than in (2.1). 
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It will be easy to prove all the inequalities stated in this section after the explicit Bellman 



functions (lL5j) and (1.6), as well as their optimizers, have been found. We provide the proofs 
in Section [81 

3. Boundary value problems and Monge-Ampere foliations 

3.1. The equations and boundary conditions. Take an interval / and split it into two 
non-intersecting subintervals: / = /_ U /-|-. For any sufficiently integrable function ip on /, 

(/M), =a-(/(V9)),_+a+(/M),^, 

where a± = |/±|/|/|. Fix two points x~, x"^ G fig such that q_x~ + a+x"*" G and consider 
a sequence of functions {v?™} such that G BMOe(/_) nBMO£(/+), {{(pn)j_^, {Pn)jj_) — 
and lim„_^oo(/('/'n))^± = Se,/(x=^). This gives 

lim {f{ipn))r = a-Bej{x~) + a+Bej{x^). 
If each (fn could be chosen so that pn £ BMO£(/), then we could conclude that 
(3.1) Bi;j{a-X~ + a+x^) > a-Bi.j{x~) + a+B^j{x'^). 

Although a_x~ + a+x^ G Og, in general 

BMO,(/) CBMOe(/-)nBMOe(/+)n{(^: {^^) ^ - {^p)] < e^} , 
which is a major difference between continuous and dyadic BMO. We will, nonetheless, enforce 



condition (3.1) for the upper Bellman function candidate B and the converse inequality for 
the lower candidate h. Thus, we look for functions B and h with the property that for all pairs 
of points x^ G such that the whole line segment [x~, x"*"] is in Q^, we have 

B{a^x'' + a+x"^) > a^B{x~) + aj^B{x'^) 

and 

b{a-X~~ + a+x^) < a„6(x^) + a-^b{x^), 

for all a± > with a+ + a_ = 1. In other words, we look for B and b that are concave and, 
respectively, convex on any convex portion of fi^. We will refer to such functions as locally 
concave and locally convex, respectively. If we also assume sufficient differentiability on B and 
b, we get differential analogs of these finite-difference inequalities: 

d^B <fb 

> 0, ^ > in 



dx^ ~ ' dx"^ 

In yet another restriction, the way Bellman function candidates are used in subsequent proofs 
suggests that we need to require that the candidates' concavity /convexity be degenerate, i.e. 
we require that 

det^^j=0, det^^j=0 inn,. 

We have natural boundary conditions for the candidates. Observe that if {p'^) j = {p)'^j, then 
p is constant on / and so {f{p))j = /(xi), giving B(xi,xf) = /(xi) and &(xi,xf) = /(xi). 



In addition, since / is assumed even, the Bellman functions (1.7) and (1.8) do not change if 
we replace p with —p in their definitions and so 

S(xi,X2) = S(|xi|,X2), 6(xi,X2) = 6(|xi|,X2). 

Accordingly, it is enough to construct candidates i3, 6 in the half-domain Q,f *== r2en{xi > 0}. 
Because of the symmetry, we impose a zero Neumann condition on the "internal" boundary 
xi = 0. 
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In what follows we will use the notation gz '= g| for any function g and variable z. Thus, 
we set out to solve the following boundary value problems for candidates B and b : 



(3.2) 



(3.3) 



BxixiBx2X2 ~ Bx^x2i Bx^xi ^ 0, Bx2X2 ^0 ill 
B{xi,xl) = f{xi), Bx^ |^j=o = 0, 

bx\X\bx2X2 ~ bx-iX2l bxixi ^ 0, bx2X2 — i 
„2n 



b{xi,Xi) = f{xi), bxi Ui=o = 0. 
3.2. Monge Ampere equations and their solutions. In this part, we state a general result 



that will help us solve equations (3.2) and (3.3). A homogeneous Monge-Ampere equation says 
that at every point in the domain the Gaussian curvature of the surface defined by the solution 
is zero, i.e. there is a direction along which the solution is linear to the second order. Such 
directions form a vector field and the integral curves of this vector field turn out to be straight 
lines. These linear trajectories foliate the whole domain and — unless the defect of the Hessian 
is allowed to exceed 1 at a particular point — do not intersect. 

There exists a formal way of obtaining Monge-Ampere solutions as functions that are lin- 
ear along certain trajectories. From this viewpoint, describing a foliation of the domain by 
such trajectories (a Bellman foliation, in our parlance) is equivalent to determining the corre- 
sponding solution uniquely. The theorem below formalizes this method of foliations. We only 
consider the case of plane domains here. The general formulation, as well as a simple proof, 
can be found in |VVlj . 

Theorem 3.1. Let Q be a plane domain and G = G{xi,X2) be a C"^ function satisfying the 
homogeneous Monge-Ampere equation in : 

i2 



Gx^X\Gx2X2 — 



X\X2 ' 



and such that either Gx^xi or Gx2X2 0- 
Let 

(3.4) ti = Gx^; t2 = Gx2; to = G - tixi - t2X2. 

Then the functions t^ are constant on each integral trajectory generated by the kernel of the 
Hessian Moreover, these integral trajectories are straight lines given by 

(3.5) xidti + X2dt2 + dtQ = 0. 



In what follows we will refer to the trajectories (3.5) as extremal trajectories, since they 
are used not only to find a suitable Bellman function candidate, but also to build extremizing 
functions/sequences proving that the candidate is, indeed, optimal. 

3.3. Main empirical principles. According to the preceding discussion, we need to look for 



a candidate G = G{xi,X2) satisfying either (3.2) or (3.3). It would be desirable to have the 



Bellman function for each / given by a single G expression that is either concave or convex 
in the whole O^. However, as we will see shortly, this is asking for too much. Instead, we 
will have our Bellman functions built out of several canonical Monge-Ampere solutions, each 
in a portion of il^, glued together so as to preserve the sign of the generalized second 
differential. A word about the nomenclature: G will stand for a generic Bellman candidate, 
whether upper or lower, whether defined on the whole domain or its part; as before, B will 
stand for a similarly generic upper candidate, and b for a lower one; when the emphasis is on 
the geometry of a canonical solution, the solution will be designated with its own name and 
necessary indices, as will be the sub-domain on which the solution is built; finally, we will omit 
the indices / and e when no ambiguity arises. 
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We will build our canonical solutions by constructing the foliations of the corresponding 
portions of the domain by extremal trajectories, while adhering to several empirical principles. 



Since the Monge-Ampere apparatus (3.4)-(3.5) does not differentiate among various solutions 



of (3.2) or (3.3), additional arguments, having to do with the nature of the Bellman function, 
are needed, chief among which is how extremal trajectories are used to construct optimizers. 
These are our principles: 

Principle 1: Symmetry. Since each Bellman function is even in xi, each Bellman foliation 
is symmetric with respect to the line xi = 0. 

Accordingly, we only need to build each Bellman foliation in the half-domain 0^, but have 
to consider carefully what happens on the internal boundary xi = 0. In which way can a 
trajectory intersect this boundary? Since the picture is symmetric, any such trajectory will 
have its counterpart from the domain n {xi < 0} hitting the same point (0, X2). When the 
domain $7^ is considered as a whole, we will have two trajectories intersecting at a point. This 
means that either they are two halves of the same horizontal trajectory or the Hessian has 
defect 2 at the point, meaning it is the identically zero matrix, and so the Bellman function is 
a linear function near that point. Indeed, we will encounter both situations below. In either 
case, since Gx^ = = when xi = 0, the function G depends only on X2 in the subdomain 
that includes the line xi = 0. 

Principle 2: Tangency. Any trajectory intersecting the upper boundary of ile, X2 = x\ + e^, 
must do so tangentially (unless the Bellman function is linear in the part of the domain 
containing the point of tangency, as in this case any straight line is an extremal trajectory). 

Let us explain this principle: observe that once the Bellman foliation for the problem is 
determined, every point of the domain can be put on one of the trajectories. That point, say 
X, prescribes the averages over the interval Q of each of the functions over which the extremum 



is taken in definitions (1.5) and (1.6). If such a function is optimal (or close to optimal), then 
when the interval is split into subintervals, Q = U the pairs of averages x^ should 
remain on the same trajectory. In such a split, x will be located between x~ and x"*". We must, 
therefore, have a tangential intersection of the trajectory with the upper boundary, otherwise, 
if we take x^ close enough to x, one of the endpoints will exit 0^. 

Principle 3. Optimality. The upper Bellman function is the smallest locally concave solution 
of the Monge-Ampere equation, while the lower Bellman function is the largest locally convex 
solution. 

This principle may be intuitively clear, and it will be rigorously demonstrated in Section |6| 
where Bellman induction on scales is used to show that any locally concave function B satis- 
fying i?(xi,xf) = /(xi) is a pointwise majorant of -Bj^g, while any locally convex solution b 
satisfying h{xi,x\) = /(xi) is a pointwise minorant of bf^^. The terms "super-solution" and 
"sub-solution," respectively, are typically used for such candidates. 

4. Local Bellman candidates: the four building blocks 

Having laid down our basic principles, we start building the foliations (and so the candidates) 
that comply with these principles. Since we know the behavior of any extremal trajectory 
touching the upper boundary, it is convenient to start with such trajectories. The tangent line 
at a point (a, + e^) is given by 

(4.1) X2 = 2axi+ - a^. 

Each such tangent intersects the lower boundary X2 = x\ at two points, 

u± = a lb e. 

Let us show that if the foliation being built includes a family of such tangents, the whole 
tangent line cannot be a single extremal trajectory, since Gx2X2 changes sign depending on 
whether x is to the left or to the right of (a, -b e^). Indeed, for every x on the tangent 
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line (4.1), ti, t2, and to are functions of a only. Recall that t2 = Gx^- Then Gx2X2 — t'2{a)a: 



X2- 



Since a is fixed along this line, t2{a) is constant. From (4.1), 

(4.2) = W7-^ r , 

2[xi — a) 

which changes sign at xi = a. 

Therefore, if the foliation contains a family of such tangents (which justifies differentiating 
with respect to a above), each extremal trajectory continues either to the right of the point 
of tangency or to the left, but not both. This is in contrast to the situation, also considered 
below, when the sub-domain being foliated lies entirely under one two-sided tangent. 

We will now consider the four sub-foliations out of which we will later build two complete 
Bellman foliations (one for the upper and one for the lower function) for various ranges of p. 
The sub-foliations are two families of one-sided tangents and two "phase transition regimes" 
used to connect those families smoothly. The four are, in order of presentation: 

(1) A collection of one-sided tangents for which the point of their intersection with the 
lower boundary is to the right of the point of tangency; 

(2) A collection of one-sided tangents whose point of intersection with the lower boundary 
is to the left of the point of tangency; 

(3) Any foliation of the convex compact set lying under a single two-sided tangent. There 
are several ways to foliate such a set, depending on its location; in each case the 
bounding tangent is an element of the foliation; 

(4) Any foliation of a curvilinear "triangle" located between two differently oriented one- 
sided tangents sharing a point on the lower boundary. As we will see, in such a triangle 
the Bellman function must be linear and so the very notion of a Bellman foliation is 
trivial, as every straight line is an extremal trajectory. 

Remark 4.1. It would, of course, be preferable to not have to combine various solutions, instead 
having a single family of tangents as the Bellman foliation on each side of the line xi = 
with a single transition regime containing that line and connecting the two families. Such a 
situation does, in fact, occur for f{s) = |s|^, p > 1. However, for < p < 1 the solutions 
corresponding to each family of tangents change concavity at various points throughout the 
domain, necessitating the introduction of other transition regimes. 



4.1. The tangents with u = a + e. Here we consider a family of tangents (4.1) with a E 
[ai,a2], or u G [ni,U2], Ui = Ui + e. The tangents under consideration are one-sided, extending 
to the right of the point of tangency, i.e. we have xi S [a,u]. The Bellman candidate built 
along these trajectories will be called (or, more fully, {x;ui,U2)) and the part of il^ 
so foliated, Qp+{ui,U2)- The foliation is shown in Figure [!} 



The candidate F+ is linear along each trajectory ( |4.1[ ) and satisfies F~^{u, u^) = f{u), thus, 
(4.3) F^{xi, 2axi + - a^) = m{u){xi - u) + f{u). 



Let us calculate t2 = F^ . Using (4.2) and the equality = o,x2-, we get 



m! iu){xi — u) — miu) + f' iu) 1 , f iu) — em! (u) — miu) 
to = ; ; = -fn iu) H z r 

Since t2 is fixed whenever a (or u) is fixed, this gives 

(4.4) t2{u) = ^m{u) 
and 

(4.5) em'{u) + m{u) = f{u). 
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Figure 1. The foliation of Qp+{ui,U2). 



Solving equation (4.5) yields 



miu) 



and so we obtain our solution, defined for x G 0^^+ (t^i, ^2), 

(4.6) F+(x;ni,n2) = e-"/^(c+- f V(s)e^/^ ds) (xi - ^x) + /(n), 
where u is given as a function of x by 

(4.7) u = n^. = xi + e — -v/e^ — X2 + xf . 



The last equation comes from plugging n = a + e into (4.1), and the minus sign in the front 
of the square root was chosen, because x\ G [a, n]. 

We will be concerned with whether (4.6 ) gives a concave or convex candidate for each specific 
choice of / (and in what parts of the domain). To that end, let us check the sign of F^^^^. 
Using (4.4) and (4.5), we obtain 

1^x2x2 = i'2{u)u^2 = \rn"{u)u^^ = ^ (/"(n) - f'{u)/e + m{u)/e) u^^. 



Since Ux2 > by (4.2), we conclude that sgn{F^^^^) = sgn(r+), where 

(4.8) T+{u) = ee''/'m"{u) = {f"{u) - /'('u)/e)e"/" + -C+\ ^ f {s)e'/' ds. 



Integrating by parts twice gives 

(4.9) r+(n) = (/" (m) - f {ui)/e)e^^'' +^C+ T (s)e^/^ ds. 



4.2. The tangents with u = a — e . We now switch to considering those one-sided tangents 
whose point of intersection with the bottom boundary curve lies to the left of the point of 
tangency. The picture is symmetric with the previous one with respect to the X2-axis. Thus, 
we replace in the preceding formulas xi by — xi, u by — n, a by —a, /(s) by f{—s), ui 
by — U2, and U2 by —ui. It is possible, of course, to derive all formulas independently, but 
we only indicate which changes are necessary. The solution constructed will be called F~ 
(alternatively, F~ {x;ui^U2)) and the portion of being foliated, VLp-{ui,U2)- The foliation 
is shown on Figure [2j 

Equations (4.3) and ( |4.4[ ) remain the same, but equation ( |4.5[ ) changes to 

(4.10) em'{u) - m{u) = -f'{u), 



SHARP LP ESTIMATES ON BMO 



11 




X2=2axi+e —a^ 

Figure 2. The foliation of Qp~{ui,U2). 

which gives 

ni{u) = e"/^ (C + - / f'{s)e-'/' ds 
This yields our canonical solution, defined for x E ^lp-{ui,U2) ■ 



F-(x;^xi,n2) = e"/^(c+- / f {s)e-'/' ds) { 

^ Ju 



xi-u) + f{u), 



with 



u = n_ = xi — e + -y /e^ — X2 + xf. 



As before, the sign of F^^^^ will be of interest. We have 

1 . 1 

-m (ujUxo = — 



Px2X2 = i'2{u)ux2 = 7:m"{u)ux2 = :^{-f"(u) - f'{u)/e + m{u)/e)ux2- 



Since Ux2 < by (4.2), we have sgn{F^^^^) = sgn(r_), where 

, . 1 1 ru2 

(4.11) r_(u) = -ee-"/^m"(u) = {/"{u) + /'(n)/e)e-"/^ " -C' " 4 / f'is)e-'/' ds. 

Integrating by parts twice gives 

1 ru2 

(4.12) r^(^u) = {f"{u2) + f'{u2)/e)e-^-/' - -C - / r{s)e-'/'ds. 

^ Ju 

4.3. The region under a two-sided tangent. Next, fix a and consider the tangent X2 = 
2axi + — 0? . It intersects the lower boundary at the points (a it e, (a it e)^) and bounds the 
convex subset of Vie 

= {xl < X2 < 2axi + — a^}. 

(In our nomenclature, each domain is designated by the name of the corresponding Bellman 
candidate; in a bit of inelegance, here we have to define and name the domain first. Naturally, 
we reserve the name L for each candidate built on ^II-) 

In order for any foliation of $7L(a) to be a part of a symmetric foliation of the whole Qe, we 
must have either a = (and so the condition Lx^ Ui=o = must come into play) or |a| — e > 0, 
i.e. nL{\a\) C nf. 

Let us first construct the solution for the set 0,^(0) = {xf < X2 < e^}. The Bellman 
candidate here will be called Lq. As explained in the formulation of the symmetry principle, 
Lq must be a function of X2 only. Thus we have ^0(2^) = 5(2^2) for some function g. On the 
boundary, Lo(xi,xf) = g{x1) = /(xi) and so g{x2) = /(y^). Thus, 

(4.13) Lo(x) = f{^). 
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The corresponding foliation of ^lii^) consists of horizontal hnes. 

We now fix an a > e and construct a solution (there are actually two) in the region flL{a). 
Assume that there is a family of trajectories foliating (a part of) this domain. Say such a 
trajectory intersects the lower boundary at the point {u^v?) and we have a whole interval of 
such values u. We do not want our trajectories to intersect the upper boundary of Q,i{a), the 
line X2 = 2axi +8^ — 0?, since we want to have that line as one of the trajectories. Thus each 
trajectory exits ^^.(a) through a point (v,w^) on the lower boundary with v = v{u). Then 
the trajectory is given by 

(4.14) X2 = {v + u)xi — vu, 
and we have, for a candidate L, 

L{xi, (v + u)xi — vu) = m{u){xi — u) + f{u). 
On the other hand, we must have L{v,v'^) = f{v), hence, 

(4.15) ™ = M^. 

V — u 

As before, let us find t2 = L^^- We have t2 = [m'{xi — u) — m + f'{u))ux2- Using (4.14) 
gives 



{v' + l)(xi — u) + u — v' 



and so 



m' m'{v — u)/{v' + 1) — m + f'{u) 
v' + 1 {v' + l){xi — u) + u — V 



which means 

(4.16) t2 = -^, m'{v-u)= {m- f'{u)){v' + 1). 

V + 1 



m! 



The last equation, together with (4.15) and a bit of algebra, gives 



f{v) - fju) , nv)+nn) , 

V = V . 



V — u I 

If / 0, then 

/(z;)-/(n) ^ r{v)^r{u) 
V — u 2 ' 

which, for f{s) = |s|P, is possible only if v = —u, i.e. only in the just-considered case of Lq. 

Therefore, v'{u) = and so the trajectories enter through different points for different values 
of u, but exit through the same point {v, v'^). Again, since we want the line X2 = 2ax\^E^ — c? 
as a trajectory, we should have the exit point in one of the corners of r2i(a), i.e. either 
V = a — E = oru = a + e = u+. In fact, each of these choices yields a solution that is either 
concave or convex in the whole region ^i^{a). 




Figure 3. The region foliated according to L 
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Figure 4. The region 0,L{a) foliated according to 

To determine whether the solution L is concave or convex for each choice of v, we compute, 
using (4.16) with v' = 0, 

r - "f \ - "^"(^) 

Thus sgn(L2^2X2) = sgn[m"(u)/(xi — i')]- A simple calculation shows that 'm"{u) = ^/"'(C) for 
some ^ between u and v. For v = U- we have xi > v, while for v = we have xi < v. 
Thus we obtain two solutions in ^l{o,) : 

(4.17) L±(x; a) = _ + /(^^), 

u± — u 

with u± and u given by 

(4.18) ii± = a±e, n = , 

Xl — ti± 

and 

(4.19) sgn(L- ,J = sgn(r ), sgn(L+ ,J = - sgn(r). 

The last two identities make sense if does not change its sign on the interval {a — e,a + e), 
as is the case for f{s) = \s\p. 

4.4. The region between two tangents. Here we fix u and consider the region Qt{u) 
between the two differently directed one-sided tangents sharing a common lower-boundary 
point (m, u^) (see Figure [5]) : 

CIt{u) = {u — e < Xl < u + e , 2uxi — + 2e\u — xi\ < X2 < x\ + e^}. 

As before, symmetry considerations dictate that for each r2j'(n) we must either have n = 
or ri-rdnl) C (that is \u\ > e). In either case, the need for a solution — let us call it 
T = T{x;u) — in this region arises when we have to glue two foliations (typically, those for 
and or L~ and L~^)- This allows us to use the optimality principle: since each of 
the two tangents bounding Qxiu) is an element of a Bellman foliation in a portion of Og, 
the compound Bellman candidate being built is linear along each tangent. Since we are to 
construct either the smallest concave or the largest convex Monge-Ampere solution in ^Ixiu), 
T is a linear function of x. 

Let us first briefly consider the case n = 0. We are looking for a function 

To{x) = aixi + 02^:^2 + ao 

on Qx{0)- As discussed earlier, Tq is a function of X2 only and so ai = 0. In addition 
ro(0,0) = /(O), which gives 

(4.20) ro(x) = ax2 + /(0). 
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The constant a is determined by reading the boundary value off the tangent X2 = 2exi; that 
value, in turn, depends on the other components of the global Bellman foliation. We will see 
in the next section how this simple step is accomplished. 

The situation when li 7^ is more involved. In theory, different choices of / may imply the 
need to use to glue various combinations of foliations of and l^^i, described in the 
earlier sections. Thus, in general we are looking for a linear candidate T in the form 



(4.21) 



such that it is equal to a particular candidate, G^, along its left bounding tangent, the line 
X2 = 2(u — e)xi — u'^ + 2ert, and to another candidate, G"*", along its right bounding tangent, 
X2 = 2(m + e)xi — — leu. That is we want to ensure that 

(4.22) T(x^;n) = G^(x^), 

where is any point on the left tangent and x~ is any point on the right one. 

The left and right tangents are assumed to be extremal trajectories for G~ and G+, respec- 
tively, and so each function is linear along the appropriate tangent. We already have equality 
at the corner (n, u^), so it is sufficient to glue our solutions at the other two corners of il7^(ii) : 



(4.23) 



G (n — e, — 2eu + 2e^) = Q!i(u — e) + a2(v^ — 2eu + 2e^) + 
G+(n + e, + 2eu + 2e^) = ai{u + e) + a2{u^ + 2eu + 2e^) + ao- 



In this paper, the only situation where we encounter r2r(n) with n 7^ 0, is when it is used 
as a transition regime between Qp+{ui,u) and Qp-{u,U2) for some numbers ui,U2. It turns 
out that such a transition places an important restriction on u. Let us elaborate. 




X2=2{u+e)xi—u —2£U 



X2=2{u—e)xi—u+2eu 



Figure 5. The region ^xiu) connecting the foliations for and F 



We are looking for a solution in the form (4.21 ) according to two requirements. The first one 
is that relations (4.23) be fulfilled, with G~ = F~ and G"^ = F~^ . The second requirement 



is that the solution must preserve the sign of the (generalized) second derivative along any 
direction. This second requirement is, in fact, the reason why T is necessary: both F^ and 
F~ are continuous in the whole domain, but their various second derivatives change signs. 
Recall that we can write F'^ and F~ as 



F^ 



m^{u){xi — u) + f{u), F = m^{u){xi — u) + f{u), 



where the coefficients ?ti+ and m_ satisfy the differential equations (4.5) and (4.10), respec- 
tively: 



(4.24) 



1 



m_ 



{m- - f). 



SHARP LP ESTIMATES ON BMO 



15 



For G- = F+ and G+ = F", dTBl) gives 



m_|_(n)(— e) + /(n) = ai{u — e) + a2(u^ — 2eu + 2e^) + oq, 
m^{u)e + f{u) = ai{u + e) + a2{u^ + 2eu + 2e^) + oq. 
Together with the boundary condition 



f{u) = aiu + a2U^ + ao 



this yields 



(4.25) 



q;2 



2e 



rri- 



4e 



772-1 



7Tl_ 



It, 



2 4e 

This determines the function T up to the parameter u, which is not free. It has to be 
chosen so that the BeUman candidate G constructed of the blocks , T, and F~ is either 
concave or convex in the combined domain Qp+{ui,u)UQTiu)^^F- ^2)- We cannot enforce 
differentiability along the tangents bounding ^It^u), but can, and will, require that the jumps 
in the first derivative(s) be of the same sign on both boundaries. We will see in Section [g] that it 
is enough to check this in a single direction transversal to the boundary. Let us do so for Gx2 ■ 
Along the left and right bounding tangents the jump in this derivative is, respectively, 02 — F^^ 
and a2 — F~^ . Recall that each and F~^ has a constant value along the corresponding line, 
given, respectively, by 



X2 



1 

Ye 



(m_|_ — /') and 



X2 



1 

2e' 



m_ 



where we have used (4.4) and (4.24). Taking into account formula (4.25) for a2, we can write 



the compatibility condition as 
m 



- m+ 1 
1^ + 2^^"^^ 



/' 



mi 



1 



f) 



1 



[m_ 



+ m_ -2/T > 0, 



4e 2e 
and, finally, as 

(4.26) m+(n) + m_(u) = 2/'(n). 

We see that, indeed, in gluing F^ and F~ via T, the parameter u cannot be chosen 
arbitrarily. In fact, only in certain circumstances such a value u exists. If it does, we can write 



the function T using (|4^, (|4^, and as 
(4.27) T{x-u) = f{u){xi-u) + 



m-{u) - f'{u) 



2e 



{-2xiu + X2 + u^) + f{u) 



This completes the description of the four basic blocks out of which we will assemble global 
Bellman candidates. The next section is devoted to this task. 



5. Global Bellman candidates 

In this section, we construct a set of global Bellman candidates, i.e. candidates that have the 
same sign of the generalized second differential in the whole domain VL^. The main emphasis 
is on our specific choice of the embedding function /, /(s) = |s|p, although some results are 
stated in more generality, which will be useful a bit later. 

Although the global Bellman foliations do significantly depend on the range of p consid- 
ered, some useful statements can be made for all p. In the previous section, we built general 
candidates F"'"(x; ui, 712) and F~(a;; tti, U2), each determined up to a constant C. To specify 
this constant for F^[x]Ui^U2), we need to know its left neighbor; for F^ {x;ui,U2), we need 
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to know its right neighbor. If has no left neighbor, i.e. if ui = — oo, the constant is 
determined using hmiting considerations, and similarly for F~ in the case U2 = oo. 

From our limited arsenal of canonical blocks, only blocks of the L type can be located 
directly to the left of F^{x;ui,U2). This gives us a boundary condition for on the line 
X2 = 2(mi — e)xi — u1 + 2uie, shared by the two canonical sub-domains. We know from (4.17) 
and ( |4.18[ ) that 

f{a-e) + f{a + e) 
2 

Setting a = ui — e and equating the result with F~^{x;ui,U2) from (4.6) for u = ui and 
xi = ui — e, we obtain 

f{ui) - fjui - 2e) 

2£ 



L{a, c? + e^; a) 



C 



which gives 

(5.1) F+{x;ui,U2) 



1 



-u/e 



f{ui) - f{ui - 2e) 

2 ^ 



(xi -u) + f{u), 



where u is given as a function of x by ( |4.7[ ). 

Let us specify this formula for /(s) = |s|^ in two cases we will need: ui = e and ui = 2e. 
In the first case we have, after a change of variable in the integral. 



(5.2) F+(x;e,n2) =peP-^e-"/^ 
and in the second one, 

(5.3) F+(x;2e,U2) =eP-^e-"/^ 
where 



u/e 

tP-^e* dt 



(xi — -u) + n^, xeQp+{£,U2), 



2P-'e^+p 



u/e 



t'P-^e^dt 



{xi—u)+uP, X £ Qp+{2£,U2), 



u = xi + e 



— X2 + x\ . 



We then have sgn(F^^2) = sgn(r+), where, according to (4.9), 
(5.4) T+{u) = p{p - 2)8^-'^ 

if ui = e, and 



ru/e 

e + {p-l) J tP-^eUt 



(5.5) 



T+{u) = {p-l)ip-2)eP-^ 



2^'-V+p 



u/e 

tP-^e' dt 



if ui = 2e. 

Now, let us consider the case when the block ^}p+ has no left neighbor, i.e. when ui = — oo. 
The argument here is subtler. Namely, we first look at the sign of F^^^^ and then invoke the 
optimality principle. For large negative u we have sgn{F^^^^) = sgn(r-i-) = sgn(C), with r+ 
given by (4.9). Thus, if C < 0, we have an upper Bellman candidate. Since we want the 



smallest upper candidate, and F~^ given by ( |4.6[ ) decreases as C grows, we set C = 0. On 
the other hand, if C > 0, we have a lower candidate, which we want to maximize; this, again, 
leads us to take C = 0. By symmetry, Clp+{—oo,U2) must be contained in ilef^i^i < 0} and 
so we can restrict ourselves to U2 < 0. Therefore, we obtain the following expression: 
(5.6) 



F+{x;- 



-00,U2) 



p^p-l^-u/e 



tP'^e'^dt 

u/e 



{u—Xi) + \uY', X £ Q,p+{ — CO^U2)^ U<U2<0, 
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with u given by 



u = xi + e — \/ e'^ — X2 + xf . 



The consideration for F (x;ni,U2) is entirely symmetrical: (5.1), (5.2), (5.3), and (5.6) 
become, respectively. 



(5.7) F {x;ui,U2) = 



u/e 



f{u2 + 2e) - f{u2) 



U2 



/ f'{s)e-''' ds 



{xi -u) + f{u), 



F (x; Ml, — e) = pe^ e 



p-l u/e 



-u/e 



and 
(5.8) 

F-(x;ni,oo) =peP-^e"/^ 
with u given by 



tP-'e' dt 



-u/e 



{u — xi) + \u\P, X G fi^- (ui, — e) 



(n — xi) + |m|^, X £ Qp-{ui, —2e), 



tP-^e-^dt 



u/e 



{xi — u) -\- uP , X G Op'- (ui, oo), < ui < u. 



ti = xi — e + a/ — X2 + xj . 



We will also need sgn{F^,_^^^{x]ui^oo)) = sgn(r_), where, according to (4.12), 

/>oo 

(5.9) r_(u) = -p{p -l){p- 2)eP-2 / tP-^e"* dt. 

J u/e 

We will now build global Bellman candidates starting "from the middle" : the only canonical 
sub-domains that can symmetrically incorporate the line xi = are riL(O) and Ot(0). We 
fix one of these and glue other canonical sub-domains to it, so as to preserve the global 
convexity/concavity of the resulting candidate. It is convenient to split further discussion 
in two parts: p > 1 and < p < 1. 

Before we proceed, let us fix the following notation for the coefficients of x\ in (5.1 ) and (5.7) 
for /(s) = \s\P : 



(5.10) m+(n;ui) = e-"*/^ 

(5.11) m_(n;u2) = e"/^ 



2£ 



^2 - (^2 + 2g)P 1 
2e 



'^^ e"^/^ + peP-^ / tP-\' dt 



ui/e 
U2/e 



tP-'e-'dt 



u/e 



< lii < ti. 



< U < U2- 



In addition, although most of our work here is with the power function, in Section [9] we will 
need the more general counterparts of (5.10) and (5.11) in the specific cases when ui = e and 
U2 = oo, respectively: 



(5.12) 
(5.13) 



m{{u) = e""/^ / f{et)e^dt 



u/e 



mL{u) = e'''' f\et)e-Ut. 

J u/e 
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5.1. The casep>l. Let us first consider the spht = (— oo, — e) U 0^(0) U r2j?+(e, oo). 



According to (4.13), the solution in ilL(O) is given by 

: Lo{x) = X 

\p{p - 2)xP/^"2 and so sgn(G^2:r 



G{x) = Lo{x) = x^/', 



hence, in that region we have Gx2X2{x) = ^PiP ~ 2)x" and so sgn(G2,2X2) = sgn(p — 2). 
In Qp+{£,oo), by ( |5.4| ) we have sgn{F^^^^) = sgn{p — 2) for p > 1. We, thus, attempt to 
check whether setting G{x) = F^{x;e,oo) in 0,p+{£,oo) will produce an acceptable Bellman 
candidate. Along the line X2 = (the shared boundary of the two sub-domains), we have Lq = 
, while (5.2) gives F'^ = for u = e. In Section pi we will verify the convexity /concavity of 



the resulting candidate in the combined domain 17^(0) ur2^+(e, oo) Subject to that verification, 
we have a complete candidate in the and hence, by symmetry, in the whole $7^ : 



(5.14) 

More explicitly, 
(5.15) M,,p(x) 
with 



M(x) 



F (x; — oo, - 
Lo(x), 
F^{x] e, oo). 



m+(u; e)(|xi| — u) + vP , 



e), x G Q.p- (—CO, — e), 
X G 17l(0), 
X G Qp+{£,oo). 

X ^VLp- (—DO, — e) U ^lp+ (e, oo) 
a; G 1^l(0), 



xi| + e - -v/e^ - X2 + x{. 



The corresponding foliation of $7^ is shown on Figure [6j This function gives an upper Bellman 




Figure 6. The Bellman foliation for the candidate M. 

candidate for 1 < p < 2 and a lower one for 2 < p < oo. Let us write separately the candidate 
for the important case p = 1, when the integrals can be evaluated explicitly: 



(5.16) 



MeAx) 



1/2 
^2 1 



e"^ — X2 + x\ \ exp ^ 



+ \Je'^ - X2 + x\ 



X G 9.p-{-oo,Q) U Vtp+{Q, oo) 
X G Ol(0). 



To get the other Bellman candidate, we consider the split = 0^+(— oo,0) U f^T(O) U 
^^-(O, oo). In f]^-(0,oo) we, naturally, set 

G{x) = F-(x;0,oo). 



SHARP LP ESTIMATES ON BMO 



19 



According to (5^), s gn(F ^^^.^) = 
the line X2 = 2exi, (5.8) gives 

F-\u=o=peP-^ 
In Ot(0), we set 



sgn(p— 2). On the shared boundary of O,t{0) and i^p- (0, oo) 



G{x) = To(x), 



where, from (4.20), 

ro(a;) = ax2. 

To preserve continuity along the line X2 = 2exi, we set a = ^eP~^T{p). Again, we postpone 
until the next section the verification that the resulting candidate is locally convex/concave in 
Or(O) U Q,p-{0,oo). By symmetry, we obtain the following global candidate: 



(5.17) 

More specifically, 
(5.18) Ns,p{x) 
where 



N{x) 



F+(x;— 00,0), X £ Clp+{— 00,0), 
To{x), X G nriO), 

F^(x;0,oo), X G Clp-{0,oo). 



m^{u; oo)(|xi| — u) + vP, x G $7^+ (— oo, 0) U $7^- (0, oo), 
|ef-2r(p)x2, xGQr(O), 



U = \Xi\ — £ + \/ — X2 + xf . 



This function gives an upper candidate for p > 2 and a lower one for 1 < p < 2. The 
corresponding Bellman foliation is shown on Figure [7| We note that this geometric description 
is accurate for all p > 1, but not for p = 1. In that case, we can again evaluate the integrals 
explicitly and thus obtain 



(5.19) 



X e Qp+{-oo,0) U il.p- {0,oo) 
■ X2, X e ^t{0)- 



Therefore, A^^^i is a piecewise linear function and the defect of its Hessian is 2 in the interior 
of each canonical sub-domain involved. Thus every straight line lying entirely in r2^-(0,oo), 
Qp+{—oo,0), or r2r(0) is an extremal trajectory for N^^i. 




Figure 7. The Bellman foliation for the candidate A^. 



20 



LEONID SLAVIN AND VASILY VASYUNIN 



5.2. The case < p < 1 . In this case, we have two Bellman candidates of a more compli- 
cated nature. As before, we build each global solution starting with either Lq{x) or Tq{x) 
(thus placing either il-iiO) or ^^^(O) at the center of O^) and then extending the solution 
appropriately to the whole Q^- 



Let us first set G{x) = Lo{x) 



in flL{0). Since sgn(Lj 



sgn{p{p - 2)) < 0, 



we are building an upper Bellman candidate. It is natural to attempt to glue a solution 
F~^{x;e,U2) to this foundation, for some U2 > 0, just as we did in the previous case. Since 
sgn{F^^^^) = sgn(r+), and, by (5.4), T+(e) < 0, this is a proper choice in that we still have 
a concave candidate. However, we cannot take U2 = oo, as we did before, since t^{u) > 
for sufficiently large u. This means that the only canonical solution we can have after all 
transient effects have dissipated (i.e. for large u, after all necessary transition regimes have 
been deployed) is F~{x; ui, oo). Again, we check the sign of F~^^^ : sgn{F~^^^) = sgn(r_) and, 
by (5.9), T_(u) < for all u. Therefore, we need a transition regime connecting the foliations 
for F^{x;e,ui) and F~{x;ui,oo). We have an obvious choice, one that was considered in 
section 4.4 T{x;^) for a specific value of ^. Condition (4.26) dictates that T{x]^) would 



appropriately glue F~^(x;e,S,) and F (x;^, oo) if and only if 
(5.20) m+(e) + m_(0 = 2pe'\ 

Slightly rewriting the integrals in (5.10) and ( |5.11[ ), and letting ^ = p,e, we can reformulate 
this condition as follows: 

r/J. roc 

dz + ei" zP-^e-^dz = 2^^"^ 
J II 



(5.21) 



zP'^e' 



Obviously, we need to have /i > 1, that is — e = — 1) > 0. We must verify that such a 
/i exists and that t+(^) < 0, meaning F^ remains an upper candidate up to the line u = ^. 
Changing the variable in (5.4) and integrating by parts twice, we rewrite this condition as 



(5.22) 



We are about to prove the existence of the solution of (5.21) satisfying (5.22). We prove 
a slightly more general result, which we will need in Secion [9} Namely, let us replace the 
condition (5.20) with 

m{_(e) + m{(e) = 2/(0, 

f \ K \ K 

where the functions are defined by (5.12) and (5.13). Rewriting the integrals, we get 



(5.23) 



\ f'{ze)e'dz + e^ / /'( 



ze)e 



In addition, recall that for a general / we have sgn{F^^^^ 
e and C = in K^, 



dz = 2f'{fie). 

= sgn(r4_), where we have set, using 



Ul 



rliu) = {f"{u) - /'(n)/e)e"/^ + \ rf'{s)e^/'ds 



f f 

After rewriting, we get sgn(r:j_(^)) = sgn(r:(_(e;u) = sgn{g{^)), where 



f{ze)e'dz-f'{iJ.e) + ef"{^ie). 



Lemma 5.1. Fix e > and let f be a thrice- differentiable function on (e,oo) satisfying 
(5.24) lim e^|/'"(se)| =oo 

s— >oo 

and either 

Case 1: f'{t) > 0, f"{t) < 0, f"'{t) > 0, Vt G (e,oo), 
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or 



Case 2: f{t) < 0, f"{t) > 0, f"'{t) < 0, Vt G (e,oo). 

Then, for each e > 0, equation ( |5.23 ) has a unique solution /x* in the interval (l,oo). Fur- 
thermore, n* satisfies 



g{fJ'*) < in Case 1, 
g{fJ'*) > in Case 2. 



(5.25) 
Proof. Letting 

poo 

h{^i)= / f\ze)e^ dz + e^f" / /'(ze)e-^ - 2e^/'(/ue), 
we se that the task is to prove that h has a unique zero n* in (l,oo). We calculate: 



h'ifi) = 2e2^ 
Thus, 



f'{ze)e-' dz - f'{fie)e-^ - ef"{i^e)e-^ 



POO 

2£2g2M / f"^ze)e-'dz. 



Case 1: h'{n) > 0, V/x > 1, and, from (5.24), lim^^oo h' (ji) = oo. In addition. 



Ml) 



/OC /"OO 
f'{ze)e~' dz - 2ef'{e) = ee^ J f"{ze)e-' dz - f'{e)e < 0. 



- Case 2: /i'(/u) < 0, V/x > 1, lim^^oo = -oo, and h{l) > 0. 

In each case, this implies the existence of a unique root fi^ of h. 
To finish the proof, observe that 

r^^ poo 

e"^* / /'(ze)e^(iz = 2/'(/x*e)-e^* / f'ize)e~'dz 

and so we have 

poo poo 

g{l2,) = -e^* / f'ize)e~' dz + /'(/x,e) + e/"(^.e) = -ee'^* / /"'(ze)e-^ dz, 



which yields (5.25). 



Setting /(s) = \s\^, we obtain an immediate 



□ 



Corollary 5.2. For each p < 1, equation (5.21) has a unique solution in the interval 
(l,oo). Furthermore, /u* satisfies (5.22). 

Remark 5.3. It is easy to show that ^*(p) — )■ oo, as p — 1^ and /x*(p) — ?■ 1, as p — )• — oo. 



From now on, let us denote the solution of (5.21) simply by /x; also let ^ = ^e. The lemma 
just proved means that we have, indeed, succeeded in building a complete Bellman candidate. 
On r2T(C) that candidate is given by T(x;^), where, according to (4.27), 



T{x; = p^-^rri + - (m_ (^; oo) - pe-^){-2x^i + X2 + C") + (1 " p)e ■ 
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Extending, as before, the solution to the left of the line xi = by symmetry, we can write 
down our global candidate: 

— cxD, — X e ^p+{—oo,—(,), 
T{x;-0, xenri-O, 



(5.26) 



P{x) 



F-(x;-e,-e), xenF-{-C,-e) 



Lo{x), 
F+{x;e,0, 
T{x;0, 
^F-{x;^,oo) 



X G Ol(0), 
X e np+{e,£,), 

X G nriO, 



This representation exhibits the geometric structure of P. In addition, we need a usable for- 
mula: 

?7i_(n_; oo)(|xi| — U-) + u^, x £ Qp+{—(X), — ^) U 17^- (^, oo), 

pe-'\xi\ + oo) - p^p-i)(-2|xiie +X2+ e) + (i - p)e^ 



(5.27) Pe,pix) 



m+{u+;e){\xi\ - n+) + u\, 



xGnri-O^^AO, 

X G -e) urij7'+(e,0, 

X G 17l(0), 



where ^ = fie, fi is the unique solution of (5.21) in (l,oo), and 

\xi\ + £ - Je"^ - X2 + xj 



u+ 



e + \ £ 



X2 + x{. 



As noted before, Pe^p gives an upper Bellman candidate for < p < 1. The foliation for this 
function is shown on Figure [Sl 




Figure 8. The Bellman foliation for the candidate P. 
To construct the lower candidate, we place ^t{^) at the center of ri^, thus setting 
(5.28) G{x) = To{x) = ax2, x G J]t(0), 

with the constant a still to be determined. We now have to glue a canonical candidate to the 
right of Tq{x). Geometrically, we need a candidate whose foliation includes the line X2 = 2£X\ 
and, thus, have three choices: -F~(x;0,/3), for some /3; L^(x;e); and L^(x;e). 
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Let us first examine F^(x;0,/3). We have sgn{F~^^^) = sgn(r_), where, according to (4.11), 
r_(u) = p{p — l)0{u^^'^) for small positive u. This means that near the boundary of O,t{0) 
(i.e. the line u = 0), we have F~^^^ < 0, that is F~ gives an upper candidate, while we are 
building a lower one. 



For the other two possibilities, (|4.19|) gives 

=F sgn(p(p 



sgn(L^ 



l)(p-2)), 

meaning only L~ gives a lower candidate. Therefore, we set 

G{x) = {x;e) = X2^^x1^^ , x G ^^^(e), 

and f{s) 



where we have used the " — " part of formula (4.17) with u 



We are now in a position to determine the constant a in (5.28): setting ro(xi,2exi) = 
L~ (xi, 2exi) gives 

a = {2ef-\ 

Having determined our candidate in ^It{0) U ili^e), we now have to glue another canonical 
solution to L~ . Observe that for sufficiently large xi we expect our candidate to be given by 
F"^(x;7,oo) (according to (5.9), its counterpart, F~ , determines an upper candidate and so 
does not work here). We attempt to take 7 = 2e, i.e. glue F^ directly to L~ without further 



transition regimes. From (5.5), sgn(r+(n)) > 0, Vn > 2e. Therefore, we have obtained the 
following complete lower candidate: 




Figure 9. The Bellman foliation for the candidate R. 



L+{x;-e), 
(5.29) R{x) = { To(x), 

L^{x; e) 
F+(x;2e,oo), 

Written explicitly, the function R is given by 

m+(n; 2e)(|xi| — u) + 



F {x; —00, —2e), x G ^7^?- (— 00, — 2e), 
X G ilL(-e), 
X G Ot(0), 
X G ^^(e), 
X G fiF+ (2e, 00). 



(5.30) 



Re,p{^) 



Xl 



|2-p 



(2e) 



p-2 



X2, 



x e ^Ip- (—00, — 2e) U Qp+{2e, 00), 
X G nL{-e)unL{e), 
X G Ot(0), 
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where 



u 



xi\+ e - \ - X2 + xf. 



The corresponding foHation is shown on Figure [9) 



6. Bellman Induction 
In this section, we first establish the local concavity /convexity properties of the global 



Bellman candidates in the whole domain 



(6.1) 



or 



(6.2) 



G{a-X + Q+x^ 



, i.e. show that each candidate G satisfies either 

> a-G{x-) +a+G{x+) 



£ such that the 



G{a-X + a+x^) < q;_G(x ) + a_|_G(x^) 

for all non-negative numbers a± such that a_ + a+ = 1 and all x^ , x^ 
entire line segment [x^, x^] is inside fig. With this in hand, we then use induction on scales to 
show that each candidate appropriately majorates (minorates) the Bellman function for which 
it was constructed. 

Before verifying (6.1 ) or (|6.2|) for our Bellman candidates, we need to make two observations. 



First, note that we can consider (6.1) or (6.2) as the statement that the derivative of G along 



every direction is decreasing (respectively, increasing) in that direction; thus, each property 
can be checked locally. We know from Section [5] that each global candidate G has the Hessian 
of the appropriate sign — and so the required monotonicity of the derivatives — in every 
canonical subdomain. Therefore, the only places where (6.1) or (6.2) needs to be verified are 
the points where the segment [x~, x"*"] intersects the boundaries between subdomains. For each 
such boundary, that verification will take the form of measuring the jump in the derivative 
along any direction transversal to the boundary. While in general our global candidates G 
are not guaranteed to be smooth at such points, we do have one-sided transversal derivatives 
everywhere and can check the sign of that jump. If that jump turns out to be at all points of 
the boundary, that means that the two solutions are, in fact, glued -smoothly. In this case, 



we will automatically obtain a global candidate satisfying either (6.1) or (6.2) in the union of 
the two subdomains. 

The second important observation is the following: whenever we glue two subdomains, it 
is always along a line that is an extremal trajectory for both local foliations. For this reason, 
while technically we need to check the sign of the jump in the derivative at all points of the 
shared boundary and in all directions transversal to it, it suffices to do so at a single boundary 
point, in a single transversal direction. Let us explain: recall that for any solution G of the 
Monge-Ampere equation, the gradient {Gxi,Gx2) = (^15^2) is constant along any extremal 
line; hence, the jump in the vector (ti, t2) is constant along the boundary and we can measure 
it at any point. On the other hand, consider new coordinates {y 1,1/2), where yi is directed 
along the shared boundary and ?/2, in a fixed transversal direction. Then Gy^ is continuous 
(so no jump in that component), because both solutions are linear functions along yi and they 
were glued continuously. Thus, we only need to measure the jump in Gy^. 

We are now in a position to prove the following 

Lemma 6.1. 



(1) 
(2) 
(3) 
(4) 



The function M^^p given by (5.15) satisfies (6.1) for 1 < p < 2 and (6.2) for p > 2. 
The function Ni^^p given by (5.18) satisfies (6.2) for 1 < p < 2 and (6.1) for p >2. 
The function P^^p given by (5.27) satisfies (6.1) for < p < 1. 
The function R^^p given by (5.30) satisfies (6.2) for < p < 1. 



Proof. According to the preceding discussion, we need to check all boundaries between sub- 
domains for each candidate. Since such boundaries are never parallel to the X2-axis, we can 
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choose X2 as our transversal direction in all cases. Therefore, to prove the lemma, we will use 
the following procedure: for each specific candidate G and each boundary between two of its 
subdomains, pick one point x = (xi,X2) on the boundary and verify that t2{x~^) < t2{x~) 
for (6.1) and t2{x~^) > t2{x^) for (6.2), where t2{x^) = Gx2{xi, X2 ^ 0). 



In all cases, by symmetry it is sufficient to check only those subdomain boundaries that are 



in nf. 



(1) In the case of M^^p, we check the jump in t2 along the line X2 = s'^ separating $7^(0) 
and Qp+{e,oo). From (|4.4[) and (|4.5[), we have, for any point x on this line. 



1 



1 



t2{x^) = ^m{e) = —{pe 



p-i 



m(e)). 



where m(e) = from (5.10) with u = ui = e. Thus, we have t2{x^ 



On the other hand, from (5.15) t2{x ) 



2^ 



2^ 



Therefore, 

the derivative jump is zero for any p and the first statement of the lemma is proved. 
(2) For Ne^p, we check the jump in t2 along the line X2 = 2exi separating ^t{^) and 
r2p-(0,oo). From (4.4) and (4.10), we have, for any point x on this line, 



t2{x~ 



1 



m'(0) = — m(0) 



2e 



where m(0) = pe^ ^ 
Thus, we have t2[x~) = 



^ dt = pe^ ^^{p) from (5.11) with u = and U2 = oo 



PpP-2 



Tip). 



On the other hand, from (5.18) t2{x 



PpP- 

2^ 



zero and the second statement is proved. 
(3) For Pe^p, we need to check three boundary lines: X2 



^r(p). Again, the derivative jump is 
between and ^lp+{e,S,); 



X2 = 2(^ + e)a;i-2^e-^2 between np+{£,^) and i^riO; and X2 = 2{S,-£)xi + 2S,e-S,'^ 
between 0,t{£,) and Qp-{^,oo). The first verification is the same as in part ([T]) above. 
The second and third are automatic: the value of ^ = /ie in (5.27) was chosen according 



to Corollary 5. 2l which ensured that condition (5.20) is satisfied. That condition, in 



turn, was a criterion for having zero jump in t2 across each bounding tangent of ^IxiO- 
(4) For Rs,p, we have two segments of the same line to check: for < xi < e, the line 
X2 = 2exi separates f^T(O) and ^^^.(e); for e < xi < 2e, the same line separates ^^(e) 
and (2e, oo). 

and t2{x 

{2eY^-p) > 0, 



In the first case, t2{x^ 
X2 = 2exi, we have t2{x~ 



i2e)P- 
: ip-l){2e)P- 



{p—l)x2 '^x'l ^ from (|5.30|). When 



and so t2{x~^)—t2{x 



which is consistent with being a lower Bellman candidate. 

In the second case, t2{x~^) is given, similarly to part ([T]), by t2{ 
^(p(2e)P-i - m{2e)), where m(2e) = {2s)P-^ from ([sloj) with u - 



X 



t2{x+) = {p- l)(2e)P-2. On the other hand, t2{x-) 
again we have smoothness of Re,p- 



{p-l){2ey 



= im'(2e) = 
= ui = 2e. Thus, 
^, as before, and 



□ 



We will make use of the following geometric result, whose proof can be found in jSVj . 



Lemma 6.2. Fix e > 0. Take any 5 > e. Then for every interval I and every ip S BMOe(/), 
there exists a splitting / = /_ U /+ such that the whole straight-line segment with the endpoints 
x^ = ^{{p) , {ip"^) is inside ^Ig. Moreover, the splitting parameter a+ = can be 

chosen uniformly {with respect to p and I) separated from and 1. 



We need another simple lemma that says that cutting a BMO function off at a given height 
does not increase its norm, which is implicitly contained in |SV| as well. 
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Lemma 6.3. Fix (p G BMO(M") and c,d gM. such that c < d. Let (fc^d be the cut-off of ip 
at heights c and d : 



(6.3) '^c,d{s) 



c, if ^{s) < c; 
(p{s), if c < ip{s) < d; 

d, if (p{s) > d. 



Then 



and, consequently, 



{'fed).; - {'Pc,d)j <{^)j- V cube J, 



IV'c.dllBMO < llf/^llBMO- 



Proof. First, let us note that it is sufficient to prove this lemma for a one-sided cut, for 
example, for c = — oo. We then get the full statement by applying this argument twice. 
Indeed, if we denote by Cdf the cut-off of f from above at height d, i.e. Cd^p = ^-oo,di then 
^c,d = -C-c{-Cdip). 

Take a cube J and let Ji = {s £ J: if{s) < d} and J2 = {s £ J : ip{s) > d}. If either Ji = 
or J2 = 0, the statement is trivial. Thus, we may assume that 7^ 0. Let = |Jfc|/|J|,A; = 
1,2. We have the following identity: 

- i^)]] - [{iCd^f)j - {Cd^fj] 

which proves the lemma, because ((/?)j^< d < {p)j^- D 

The following is the main result of this section. Its statement is similar to — if much more 
general than — that of Lemma 2c from |SVj . The proof, by induction on pseudo-dyadic scales, 
is somewhat streamlined compared to that in [SV], although its main ingredients are the same. 



Lemma 6.4. Fix e > and let B and b be two functions defined and continuous on 



for some 6 > e. Assume that B has property (6.1) and b has property (6.2) on 0.^. Let 
W{t) = B{t,t^), w{t) = b{t,t^). Lf either W or w is unbounded at —00 or +00, assume it is 
monotone for t sufficiently close to —00 or +00, respectively. 

Fix a point x E and an interval Q and take any function if E BMOe{Q) such that 
{{^)q,{^^)q) = x. Then 

B{x) > {W{^))q, 
bix) < {wi^))^, 

including the possible infinite values on either side of each inequality. 

Proof. We will only prove the statement of the lemma concerning B, as the part concerning b 
is virtually identical. We first establish the result for those (p G BMO^IQ) that are bounded 
and then approximate arbitrary BMO functions by appropriately chosen cut-offs. 

Take ip G BMOe(Q) n Observe that p £ BMO^iL) for any subinterval / of Q. We 

now build a binary tree D{Q) of subintervals of Q, where every interval / G D{Q) is split 



into two subintervals L± G D{Q) according to the rule from Lemma 6.2 The set of intervals 
of the n-th generation will be denoted by Dn{Q), so Dq{Q) = {Q}, Di{Q) = {Q±}, etc. For 
every interval I G D{Q), let x^ G Q.£ be the corresponding Bellman point, x^ = (v^)/) • 

Let x^^^ denote the step function from Q into r^^, defined by the rule x^'^\t) = if t G /, 



/ G Dn{Q) ■ Since Lemma 6.2 provides for the value of a+ uniformly separated from and 1 
on every step, we have 

max {|/|| — )• as n — )• 00. 

/ez?„(Q) 
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By the Lebesgue differentiation theorem, we have x^'^^t) — )• {(p{t), (p'^{t)) almost everywhere. 
Since (p is assumed bounded, {x^")} is a sequ ence of bounded functions. 



For each of the splits prescribed by Lemma 6.2, the line segment connecting x^~ , x^ , and 
lies in Using the property ( |6.1[ ) of B repeatedly, we get, for any n > 1, 

\Q\Bix^) > \Q+\B{x'^+) + \Q.\B{x^-) 



> \I\B{x^)= I B{x'^''\t))dt. 

/eD„(Q) •'^ 



Since B is continuous on 17^ (and thus on Oe), the dominated convergence theorem applies 
and taking the limit as n — )• 00 proves the lemma for bounded (p. 

Take now an arbitrary ip G BMO(r(Q). For c, d G M such that c < d, let (pc,d be defined 
by (6.3). We have ipc,d G L°°{Q) and, by Lemma 6.3, ||</?c,(i||BMO ^ £• Therefore, 

B{{Pc4)Q,{^ld)Q)>{W{iPc4))Q- 

We now take the limit in this inequality as c — t- —00. Since B is continuous, the limit of the 

left-hand side is B{{ip_oo,d)Q,{v>-oo,d)Q)^ where ip^oo,d '= dx{^>d} + ^X{>f<d}- On the other 
hand, by the monotone convergence theorem, the limit of the right-hand side is {W{ip-oo,d)) q- 
The same argument works if we let d — )• 00, which completes the proof. □ 

As an immediate corollary, we obtain the following 

Theorem 6.5. For any x £ ^If,, we have 

Forp>2: N^^pix) > Be,p{x), Me,p{x) < be,p{x), 

Forl<p<2: M,^p{x) > B,^p{x), iV,,p(x) < 

For < p < 1 : Pe,p{x) > Be,p{x), Re,p{^) < be,p{x). 

Proof. Let Bf, stand for any of the upper candidates in the statement of the theorem and 6^ 
for any of the lower candidates. Since each B^ and each is continuous in e, it is sufficient 
to prove that 

Bsix) > Be,p{x), bs{x) < be^p{x), Vx G ilg, 

for all 6 > £, and then let 5 — )■ e. 

Observe that each Bs and each bs satisfies the conditions of Lemma 6.4 and Bj{t,t^) = 
bs{t,t^) = \t\P. Therefore, 

Bs{x)>{\pnQ, bs{x)<{\pnQ, 

for all if G BM0e(<5) with {{(p)^, {pP') = x. Now, take the supremum over all such p in the 
first inequality and infimum in the second. □ 

7. Optimizers and converse inequalities 

In this section we construct, for each Bellman candidate G built in Section [5] and each point 
X G a test function ip^ on (0, 1) with the following three properties: 

(a) ('/'x>(o^^)= xi, {pD ^^^^^= X2; 

(7.1) {h) 99, GBMO,((0,1)); 

(c) (/(v^.))(o^,)=G(x). 

We will call each such function px an optimizer for G{x). We will often need optimizers for 
points on the top boundary of U.^, the parabola X2 = x\ + . It is convenient to parametrize 
these by the horizontal coordinate: let us denote p(^x]_,x\+e'^) i^xi ■ For the rest of this section. 
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we will simply write BMO^ for BMO£((0, 1)). If we need BMO over another interval, we will 
write the interval explicitly. 



If the candidate G corresponds to the Bellman fmiction G, the existence of ip satisfying (7.1 ) 
would immediately imply that 

G{x) < G{x), Vx G if G is an upper Bellman function, 
G{x) > G{x), Vx G Qe, if G is a lower Bellman function. 



In general, whether an extremal function exists for a particular inequality is a deep question, 
which sometimes is more difficult to answer than to prove the inequality itself. However, if the 
Bellman function for the inequality is known, one has a definitive answer as to the existence 
and the nature of optimizers. If one exists, it is found using the geometry of the Bellman 
foliation, as done in this section. If one does not exist (which is the case, for example, for all 
known Bellman functions for the maximal operator), optimizing sequences are built along the 
trajectories of the foliation. 

Since our global Bellman candidates are built out of four canonical local blocks, F~^, , L, 
and/or T, each with its own geometry, it is natural to attempt to construct a corresponding set 
of canonical optimizers, one for each local candidate. If one has such a set, one can demonstrate 
that (the appropriate line of) (7.2) holds for the global candidate G simply by showing it for 
each canonical sub-domain of il^. 

Thus, we fix a canonical block D and attempt to build an optimizer ip^ for each point 
of the domain Q^. Expectedly, ipx will, generally, depend not only on x, but also on the 
placement of Qd within il^. This is so because our local candidates, for the most part, them- 
selves explicitly depend on that placement. Somewhat more subtly, the optimizers turn out 
to depend, in some cases, on exactly how the local candidates are glued together to produce 
a global one. In Section |4j we built blocks of two kinds. Some were completely determined 
by the parameters of their domain, namely F~^{x; —oo,U2), F~{x;ui,oo), and all blocks of 
the L type; let us call these blocks complete. Others, namely F+(x;ui,U2) for ui > — oo, 
F~(x;ui,U2) for U2 < oo, and T, had undetermined constants that were only found in Sec- 
tion [5| where the neighbors of these blocks in the global context were examined; let us call 
such blocks incomplete. Furthermore, different incomplete blocks may or may not require the 
knowledge of both neighbors. Thus F~^ is left-incomplete, unless ui = — oo, since it requires 
the knowledge if its left neighbor (which, from among our limited supply, can only be an L 
block); is right-incomplete, unless U2 = oo, as it requires the neighbor on the right (again, 
an L block); and T is both left- and right-incomplete, as it requires both neighbors (these can 
be various combinations of L and F blocks). 

This division has exact parallels in this section: to determine an optimizer for a complete 
block, we only need to know its domain. For example, it is meaningful to say that a given 
function is an optimizer for L~^{x;a). However, to determine an optimizer for an incomplete 
block, we need to know its neighbors, and their optimizers, on the left and/or on the right, all 
the way to a complete block. Thus, we cannot say that a given function ipx is an optimizer 
for F^{x;ui,U2), because F^ is incomplete. We could, however, say that a function is an 
optimizer for the sequence L{x;ui—e) — )• F"'"(x; ui, U2), because such a sequence determines 
the candidate completely. 

In such situations, we will talk about optimizers for a block (alternatively, local candidate) 
in the context of a sequence. Rather than defining this notion in general, let us list the specific 
sequences we will encounter, including those consisting of a single block: L, F~^{x; —oo,U2), 
F-{x;ui,oo), L F+, F' L, L ^ T ^ L, F+ ^ T ^ F~ , and L F+ ^ T ^ 
F^ (x;ui,oo). In each case, all arrows point to the block for which an optimizer is being 
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sought. Expectedly, each such sequence starts with a block that is complete on the left and 
ends with one complete on the right. 

The main empirical consideration in building optimizers for Bellman foliations is that for 
each point x the whole construction of the optimizer ip^ should be taking place along the ex- 
tremal trajectory passing through x. This means that ipx should be such that when the interval 
/ = (0, 1) is split into two subintervals, / = /_Ul+, the Bellman points = {{cpx) 
are also on the trajectory. Indeed, for the true Bellman function, each inequality in the Bell- 
man induction of Section [6] is, in fact, an equality; thus, each split must be happening along 
the trajectory where the candidate is linear. Conversely, if m, n are two points on such a 
trajectory and we know the optimizers ipm and ipn, then the optimizer for any point x on the 
trajectory that is between m and n can be obtained simply by concatenating the two known 
optimizers: 



t 

7 

1-7 



t G (0,7), 
i G (7,1), 



where 7 = (xi — ni)/(mi — ni). 



Remark 7.1. In general, concatenating two functions from BMO^ is not guaranteed to produce 
another BMO^ function. One or both of the functions being concatenated, ipn and ipm in the 
formula above, may have to be rearranged to ensure the smallest possible BMO norm of the 
resulting function. 

Remark 7.2. As we will see, the optimizers built in this section do not depend on the choice of 
the boundary function /, but only on the geometry of the canonical subdomain and, in some 
cases, on how such subdomains are glued together. 

7.1. Optimizers for Lq{x) . The foliation of the domain ^2^(0) that corresponds to the Bell- 
man candidate Lq given by (4.13) consists of horizontal lines. For each point x £ ilii^), the 
extremal trajectory through x intersects the boundary of fig in two points, {—^/x2, X2) and 
(y^, X2). Since the only test functions available on the boundary of il^ are constants, we 
already know the optimizers for these points: ,j,2)(t) = ±y^. Therefore, to construct 

the optimizer for the point x, we concatenate the two boundary values in the appropriate 
proportion: 



(7.3) 

where 
(7.4) 



a 



X2, t G (0, a), 
x^, t £ (a, 1), 



We now verify (7.1) for ip 



Lemma 7.3. The function (fx given by (7.3) and (7.4) is an optimizer for Lq{x). 



Proof. For part (a) of (7.1), we have {ipx 



'(0,1)" 



-a^fx2 + (1 — OL)^fx2 = x\ and (y?^) 



(0,1)" 



aX2 + (1 - (y)X2 = X2- 

To show (b), we need to show that for any subinterval / of (0,1) the Bellman point x^, 
corresponding to / and p)x , is in 0^ . Observe that x^ is a convex combination of two boundary 
points, (±y^, X2). Thus, it lies on the line segment connecting these points and, furthermore, 
belongs to the convex subset r2L(0) of fi^. 

For (c), we trivially have (/(¥'xO)(o,i)= fiV^) = Lo{x). □ 
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7.2. Optimizers for 



x: a 



a 7^ 0. Recall that for each choice of a the candidates 

C( 

^ The Bellman foliation 



(re; a), given by ( |4.17| ) and ( |4.18[ ), are built on the domain ^^(a) consisting of the (closed) 
portion of lying under a two-sided tangent X2 = 2axi + — a 
for is the collection of straight lines connecting the corner point {u+,u'^) to points {u^u^) 
with u„ < n < while the foliation for consists of lines connecting to 
For each x G VIl{o), except the two corners, there is a unique such u, given by ( |4.18| ): 

X2 - vxi 

(7.5 u= , 

Xl — V 

where v stands for either u+ or u-, as appropriate. 

Again, we already know the optimizers for the points {v, v^) and (u, v?) : (p(^^^^2-j{t) = v and 
(/9(„ „2)(t) = u, respectively. Therefore, to construct the optimizer for the point x lying on the 
line connecting v and u, we concatenate these two optimizers: 



(7.6) 

where 
(7.7) 



a 



Xl — V 



u 



V 



I, t £ (0, a); 
t G (q, 1), 

{xi - vf 
X2 — 2vxi + v"^ 



We now verify (7.1) for cp 



Lemma 7.4. The function ip^ given by (7.5), (7.6), and (7.7) is an optimizer for L (x;a), 
ifv = u^, and for (x; a), if v = 



Proof. For (a), we have (992 



au + {1 — a)v = xi and 



au'^ + (1 — ct)v'^ 



X2- 



^'(0,1) ' "■^'^ ^^^'(0,1) 

For (b), we proceed as before. Let x^ be the Bellman point corresponding to tpx and a 
subinterval / of (0,1). This point is a convex combination of the points {v,v'^) and (u,ti^), 
which lie in the convex subset ^^(a) of O^. Therefore, x^ belongs to Og. 
Finally, we check the optimality relation 



(/(^)> 



(0,1) 



af{u) + {l-a)f{v) 



) of ([TTj): 



u 



(xi -v) + f{v) = L{x;a), 



according to (4.17). Here L stands for either or L , as appropriate. 



□ 



7.3. Optimizers for L(x; ui—e) 



F+(x;mi,M2) 



-u/e 



f{ni) 



F^{x; ui,U2) • Recall formula (5.1 ) for F'^ in this context: 



f'{s)e'/'ds 



Ml 



(xi -u) + f{u), 



with u given by (4.7): 



u 



xi + e 



- X2+ Xl- 

The extremal trajectories for are one-sided tangents to the upper boundary. According 
to the discussion in the introduction to this section, to be able to construct an optimizer for 
any x G ^p+{ui,U2), we first construct one for each point of the upper boundary of this 
sub-domain, {{a,a^ + £^)) cl G (01,02]}, ai = Ui — e. Then, we concatenate the optimizer on 
the upper boundary with the constant optimizer on the bottom boundary via the extremal 
tangent through x. 



Fix an a G (01,02] and consider the tangent X2 = 2axi -|- e 



intersecting the lower 



boundary at the point (n,ii^). Since all our extremal tangents are one-sided, we do not seem 
to have another optimizer with which to concatenate the constant function u. However, we 
circumvent this difficulty with the following approximating procedure. Fix a small number 
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A and let a 



e-A 
e+A 



' ^ ~ e+A' ■'^o'^ consider the point (a — A, + — 2aA), also on the 
tangent. Assume for a moment that we know the optimizer p for this point and set 

(7.8) '^("■«- 

[u, tG{f3,l). 

To get the optimizer p, we draw a tangent through the point (a — A,a^ + — 2aA), which 
intersects the upper boundary at the point (a — 2A, (o — 2A)^ + e^) and the lower boundary, 
at the point {u - 2A, (u - 2A)2) (see Fig ure 10^. If we knew the optimizer ipa— 2 A, we could 
again concatenate the optimizers on the two boundaries: 

'V'a-2A(^t), tG(0,l-A/e), 
2A, tE(l-A/e,l). 



(7. 



Pit) 



u 



Combining (|7.8|) and (|7.9|), we obtain a recursive approximation for "ijja '■ 

>a-2A(^), te(0,a), 

Mt) - lu-2A, 



u 



t G (a,/3), 




Figure 10. The construction of 
Repeating this procedure k times, we get 



(7.10) 



def 



M-2A;A, t G (a'=,a'=-i/3), 

M-2(A;-1)A, te{a^-^p,a^~^p), 



2A, 



f G (a/3,/3), 



We seem to have a major problem: we do not know V'a-2fcA- We do, however, know ijja^. 



Setting A = "^fc^ = ^ 2k^ ' obtain a known function on the top line of (7.10). Now it is 



time to let A; — t- oo (i.e. A — )■ 0). To determine ipa = \mij^^^ we write down a simple 
differential equation: take a j, 1 < j < A;, and let t = (3 be a generic point in (0, 1). Then 

Mt) - M»t) ^ 4'\t) - i,^^\at) « (n - 2iA) - (n - 2(j + 1)A) = 2A. 

On the other hand, 



V'a(t) - Mat) « Mt)t{l - a) « <(t)t 



2A 



Combining the two approximate equalities and solving the differential equation, we obtain 

Mt) = D + e\ogt. 
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Observe that 



lim a = lim ( 1 



2A 

e + A 



(a-ai)/(2A) 



,(ai-a)/£ 



Altogether, (7.10) becomes 
(7.11) 



>ai [e^^-^'^^H] , t G (0, e(''i-'*)/") , 
D + e\ogt, t E (e(''l-'^)/^ 1) . 

How do we determine the constant D? Recall that we a priori have {il^ai) ^•^= o-i and must 
ensure ii^a) -^^•j= This gives 

^^g{ai-a)/e + _ ^) _ g(ai-a)/£\ ^ e^a^-a)/e^^ - ffli) = O 



or 



D = a + e = n . 



Having constructed an optimizer for each point on the upper boundary, we are in a position 
to construct one for any point x E Qp+{ui,U2)- As planned, we consider the extremal tangent 
through X, and concatenate the optimizer ipa for the upper boundary and the constant u for 
the lower boundary. Specifically, we have 



Using (7.11), we obtain the complete expression for the optimizer: 
(7.12) 



fV'ai(7i;), tG(0,H 



u + elog^, te{iJiy,fi), 
u, t E (/X, 1) , 



where 
(7.13) 



u — Xl 
e 



, = e*^"^ and u = u^= xi + e — \/e'^ — X2 + x^ 



Let us now recall the expression for ^ai, the optimizer for L{ai,a1 + e^;ai). It is given by 
either (7.3) or (7.6), in each case with a = 1/2 : 



V'ai(t) 



ai — e 
ai + e 



ui-2e, tE (0,1/2), 
ui, te (1/2,1). 



Therefore, formula (7.12) can be rewritten as 



(7.14) 



'ui-2e, tE (0,^1^/2), 

ui, t E (/xi^/2, nu) 

n + elog^, te {fJ-iy,f^), 

u, tE(/x, 1). 



We can now prove 



Lemma 7.5. The function ip^ given by (7.13) and (7.14) is an optimizer for 
L{x; ui —e) — > F+(x; ui,U2)- 



SHARP LP ESTIMATES ON BMO 



33 



Proof. First observe that 
t 

log — dt = log u — fi + ^u, 



Now, using (7.13) and the fact that {i^ai 



log^ -dt = - 

)(o,i)="i'^^«i 



-fjLU log^ V + log V + 2(/i — ^v). 



(0,1)" 



a\ + we get 



(0,1) 



— (ui - 2e) + ^■ui + u{X - l^i'j + e(-/iz^logz/ - /" + /"i^j 



^v{ui — e) + n(l — /ui^) — ^u{ui 

Xl 



xi) + e/ii/ 



and 



<^'){o,i) 



■ X2. 



-^v log^ V + log u + 2iJL — 2^v) 



Proving that (fx £ BMOe is more delicate. Let q{t) = u + e log{t / fi) . Observe that q E 

*== logt and any interval {c,d) C (0, oo) a direct calculation 



BMOe((0,oo)). Indeed, for l{t) 
yields 



(7.15) 



'{c,d) 



a? 



1 



cd 



id- 



< 1, 



and so I G BMOi((0, oo)), immediately implying the result for q. Moreover, setting c = 
in (7.15) shows that the Bellman point corresponding to q and any interval of the form (0,7) 
is always on the upper boundary of fig. The last bit of information we will need about q, one 
that can be shown by a computation similar to that at the beginning of the proof, is that 



(7.16) 



(q) 



(0,^1/) 



ai. 



2 , 2 
aj + e . 



Next, since ipx is the cut-off at height u of the function 

'ui-2£, t£{0,iJ,u/2) 



Vx{t) 



def 



til, 

Qit), 



t G {fiu,l), 



it suffices to show that r]^ G BMO^, as Lemma 6.3 will then imply the result for (fx- Thus, we 
aim to show that (5c,d *== (f?^^ _ \2 ^ ^2 



(c,d) 



{rix)^^^^ <e' for all (c,d) C (0,1). 



We note that the only subintervals (c, d) that need to be considered in detail are those with 
c G (0,/uz//2) and d G (/uz/, 1). Indeed, if < c < d < jiv, then < e^, since ipai S BMO^, 
as has already been shown. If ijlu/2 < c < d < 1, then, again, < e^, since r]x\[^v/2,i) is the 
cut-off, at height ui, of the BM0£((/Uz//2, 1)) function q and so Lemma 6.3 again applies. 
Thus we focus on the case < c < /iz^/2 < iiv < d < 1. Let z~ = ((^x)^q^^, (^x)(oc)^' ^ ~ 



be the three Bellman points corresponding 

to the intervals (0,c), (0, d), and (c, d), respectively. Since f?x|(o,c) is constant, is on the 
lower boundary of Vte- To locate z, we note that by (7.16) we have 

and, therefore, 



for all d > fiv. As noted above, z must be on the upper boundary of rig 
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Consider now the line through z~ and z. Since z is, by construction, to the right of the 
point (ai,a^ + e^), this hne hes above the top boundary of ^^{ax) and so first exits f^e and 
then re-enters it at z. Since z is a convex combination of z~ and z"*", z"*" is located on the 



same line, to the right of 2, and, therefore, inside ilg (see Figure 11). 




Figure 11. The mutual location of 2: , z, and 



For part (c) of (7.1), we have 



3(^-")/7(s)d5 + /(^x)(l-^) 



1 [/(^xi - 2e) - /(m)] - ;u re(^-")/V'(s) + /(n) 



Ul 



F+(a;;ni,n2). 



Here we have used the formula (7.14) for on the first step, the change of variables s = 
u + elog(t/;u) on the second, integration by parts on the third, and expressions (7.13) for 
and v on the third and fourth. This completes the proof of the lemma. □ 



7.4. Optimizers for F^{x; —oo,U2). This case follows directly from the previous one. Here 
the candidate is given by 



F^{x; —00, U2) = -e 



-u/e 



(xi -u) + f{u), 



with u given, as before, by 



u = u-^- = xi + e — \ / e'^ — X2 + xf. 



We readily obtain an optimizer for this case by setting u = in (7.14): 
(7.17) 



li + elog^, t£{0,n), 



u, 



te (^,1) 



Lemma 17.51 can be restated for this case: 



Lemma 7.6. The function ip^ given by (7.17) is an optimizer for F~^{x; —oo,U2)- 

7.5. Optimizers for sequences containing F^ . To construct optimizers for F^ we use 
the symmetry F~ {xi,X2]Ui,U2) = F+(— xi, X2; — U2, — ui) . Therefore, in the optimizers given 
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by (7.14) and (7.17) we have to replace ip by 

'u2 + 2£, 

(7.18) ^,it) = 



U2, 



elog 



-If and ui by —U2. This yields 

t G {0,fiu/2) , 
t G (/xi//2,^z^) , 
^, t G {^J.ly,^I), 
t G (/x,l), 



for F (x; ui,U2) 
(7.19) 

for F~{x] Ml, oo) 



L(x; + e) and 
V'z(i) = 



elog 



t G (0,^), 
t G (;U,1), 



(7.20) 



Here 

e 



^{u-U2)/e 



and 



e + \ - X2 + XI . 



Lemma |7.5| and Lemma |7.6| can be reformulated, respectively, as 

Lemma 7.7. 

F-{x;ui,U2) 

Lemma 7.8. 



The function ip^ 
i — L(x; U2 + e). 



given by (7.18) and (7.20) is an optimizer for 



The function ipx given by (7.19) and (7.20) is an optimizer for F {x;ui,oo 



7.6. Optimizers for sequences containing T(x; u) . Constructing optimizers for a T-type 
candidate is straightforward, since T{x;u) is a linear function in ^^(ti). Here is our strategy: 
for each point x G i^xiu), draw any straight line that intersects both straight-line sides of 
(u) , but not the upper boundary of Qe , and then concatenate the optimizers for the two 
points of intersection in the appropriate proportion. Since T is linear, the resulting function 
will be an optimizer for it. This requires knowing optimizers along both bounding tangents, 
which is consistent with the fact that T is both left- and right-incomplete. Thus, we will 
eventually need to examine the specific sequences in which T shows up in our global Bellman 
candidates, in order to write down an explicit optimizer for each case. However, most of the 



construction, as well as the verification of parts (a) and (c) of (7.1 ), can be carried out without 
specifying the left and right neighbors of T. 

Suppose r(x; u) is glued to a Bellman candidate G~{x) along its left bounding tangent and 
a candidate along the right one. Let us assume that we know optimizers for along 
their respective bounding tangents. 

For any x G i^xiu) there is always a way to draw a line through x intersecting the left 
and right tangents at the points x~ and x~^ , respectively, and such that the whole segment 
is in Qxiu). (For example, at least one of the tangents to the upper boundary of 
that pass through x will satisfy this requirement. In fact, for each point of i}T{u) other 
than the three corners both tangents will satisfy it.) Then, we can concatenate the two known 
optimizers, (p^- and iPx+, to obtain a test function ip^ for the point x : 



(7.21) 



where 



(7.22) 



<Px{t) 



t G (0,/3), 
t G (/3,1), 



Xi 



X2 



We can verify parts (a) and (c) of (7.1) directly from (7.21) and (7.22). For the averages, we 
have 

-2\ 



X2. 
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To check the optimahty of ifx, we write T = aixi + 02X2 + ckq and calculate 



(0,1) 



G-(x-)/3 + G+(x+)(l -/3) 

T{x-;u)l3 + T{x-^-u){l- 

{aix^ + 023^2 + ao)/3 + [a.ix'l + a2X^ + ao)(l 

aixi + a2X2 + ao = ^"(2;; ?^)- 



Here we have used: on the first step, the optimality of (p^- and ipx+ for and G^, re- 
spectively; on the second and third steps, the boundary conditions (4.22) for T along the two 
bounding tangents; and on the last step, the fact that x = + (1 — /3)x"'". 

Before specifying ipx for each sequence involving T{x;u), let us rewrite (7.21) in the form 
that is independent of x~ and x~^. Since each bounding tangent is assumed to be an extremal 
trajectory for the corresponding Bellman candidate (the left tangent for G~ , the right one 
for G"*"), the candidate is linear along the tangent. Thus we can write fx-, an optimizer for 
G'{x^), as a concatenation of optimizers for G~{u — e,{u — e)^ + e^) and G~{u^v?) and 
similarly for {px+ '■ 

(7.23) ifx-{t) 
where 



t G (0,a_) 
t E (a_, 1) 



U, 



t e (0,a+), 
t G (q+,1). 



u 



a+ = 1 



Observe that in the expression for cpx- we put the constant value u on the right part of (0, 1), 
while in the expression for 99^+ this constant value is on the left. This is done so as to minimize 



the BMO norm of the the resulting function ipx given by (7.21) (see Remark 7.1). Using (7.21) 



in conjunction with (7.23), we get 
(7.24) fx{t) 



n, 



, , t-/3-(l-/3)a+ 
'Pu+e\ (l-i3){l-a+) 



t G (0,/3a_), 

t e (/3a„,/3 + (l-/3)a+), 
t G (/3 + (l-/3)a+,l). 



Since x is a unique - 

(m — e, — 2eu + 2e^ 



and independent of x and x"*" — convex combination of the points 
{u, V?), and (u + e, + 2eii + 2e^), we expect the weights /3a_ and 



(1 — /3)(1 — a+) in (7.24) not to depend on how the points x were chosen. Indeed, after a bit 



of algebra, we can rewrite (7.24) as 



(7.25) 

where 
(7.26) 



X2 — 2uxi + u 

4^2 



u, 

Xl 



2e 



t G (0,/i_), 
t G (/i_, 1 - 
t G (1 -;U+,1), 

X2 — 2nxi + u'^ 



4e^ 



+ 



Xl 



2e 



We now need to check that (fx constructed according to (7.25) will be in BMOg for each 



global context of T. Among our four canonical blocks, a T block can only have or L 
blocks glued to its left and only F~ or, again, L blocks to its right. Since we have already 
constructed optimizers for all and L blocks, we could proceed in generality and prove 
that an optimizer of the form (7.25) could always be rearranged — separately on (0, /i_) and 
(1 — 1) — so that the resulting function is in BMOg. However, we choose here to be more 
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explicit and consider specific optimizers for the specific sequences in wliich T appears. We 
have three such sequences an so split further presentation in three parts. 



7.6.1. Optimizers for F'^{x;— CO, 0) — >To{x)< — F (x;0,oo). This sequence appears in (5.17). 



For this case (7.25) gives 



(7.27) 



elog 



0, 



-elog 



t G - 
t G (l-/i+,l). 



To show that (fx G BMO^, take an interval (c, d) C (0, 1). First, observe that if c G [/U_, 1— ^+], 
the Bellman point x^^''^^ is the same as the one for the cut-off at height of the function 
elog(t//i_), which has been shown to be in BMOg; therefore, x^'^''^^ is in fig. The same 
reasoning applies when d G [/U-,1 — /U+J. Therefore, we will assume that c G (0,/i_) and 
dG (l-/x+,l). 

Recall representation ( 7.22[ ), which gives x as a convex combination of x^ and x^. In our 
standard notation, x~ = x^^) and x~^ = x^^'^\ We first would like to determine the location 
of the point x^^'^\ We know that x^^'^^ is the Bellman point for the interval (0,c) and the 
logarithm elog(t//i_) and we have already seen that every such point is on the parabola 



X2 



x^ + e^. Therefore, above the line m through x and x 



(see Figure 12). Since 




(0,0) 



Figure 12. The location of the point x^"^'^) 
is a convex combination of x^^'^^ and x^^'^\ we conclude that below m. Similarly, 



the point x^^''^^ (not shown in the figure) and, therefore, the whole line segment [x^'^'^^ 
is below m and so in Q^. The point x^^''^^ is a convex combination of x^"^'^^ and x'^^''^\ which 
means that it is on the line segment and, thus, in fig. 



7.6.2. Optimizers for Lq{x) — > F"'"(x;e,^) — > T(x;^) < — F (x;^,oo). This sequence ap- 
pears in (5.26). Note that we have to include in the sequence the block Lq to the left of F^ , 
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since is left-incomplete. For this case (7.25) yields 

'-£, t G (0,/i_z./2), 

C + e log 



(7.28) 



t G (/i_z//2, /i_z^), 

t G 1 - 
te (l-^+,l), 



where = and are given by (7.26). 

Instead of (fx, it is convenient to consider ifx 



elog 



-elog 



t G (0,^_z./2), 
t G (/x_z^/2, /i_z^), 
t G {^^u,fi.), 
t G 1 - /i+), 



Clearly, ipx G BMOg <^=^ c^^: G BMO^. Let us take a subinterval {c,d) of (0,1). Then 
'fx\{i_i_iy,i) is just (the appropriate restriction of) the optimizer from the previous subsection, 
given by (7.27), and so in BMOg. Furthermore, ^x\{^_u/2.i) is the cut-off of that optimizer 
at height e — ^, restricted to (//_z^/2, 1), and. hence, it is also BMOg. Therefore, in proving 
that the Bellman point x^'^'"'^ for (fx is in we only need to consider < c < ji-v /2. Where 
should we place (P. If d < 1 — /i^-, x^'^''^^ can be seen to be a Bellman point for the optimizer 
for the sequence L — )• F"*", given by (7.14) with ui = + e and n = 0; therefore, it is in fig. 
Thus, from now on we will assume that c < ii-v/2 < 1 — < d. 

As before, we appeal to representation (7.21)-(7.22). The only difference now is that the 
point x^"''^) = (— ^ — e, (^ + e)^) will not be on the upper parabola X2 = + e^, but instead on 
the lower parabola X2 = x\. Since ^ > e, this point is above the line X2 = —2exi. Therefore, 
it is above the line through x~ and x^ (let us again call it m). Since x~ = x^'''^^ is a 
convex combination of x^^'^^ and x^'^'^\ we conclude that x^'^'^'^ is below m. From here the 
consideration is identical to the one in the previous subsection and we conclude that x^^''^^ G 0,^. 

7.6.3. Optimizers for L^{x;—e) ■ 



From (7.25) we obtain 
(7.29) 



To{x) 



-2e, 



0, 
2e, 



< — L {x;e). This sequence appears in (5.29). 

t G (0,^_/2), 

t G (^-/2,l-^+/2), 

tG (l-^+/2,l). 



To prove that ipx G BMOe, take an interval (c, d) C (0, 1). If c > )U_/2 and/or d <1 — //+/2, 
the Bellman point x^^''^^ is a convex combination of one of (±2e,4e2) and (0,0). The point 
rf.(c,d) ^j-^en ig on one of the tangent lines X2 = ±2exi and, thus, in $7^. Let us assume that 
c < /i-/2 and 1 — /i_|_/2 < d. More fully, let us write c < fi- /2 < < 1 — < d, where 

/3 is given by (7.22). We have x~ = x^^'^^ as a convex combination of x^'^'^^ = (— 2e,4e^) and 
Therefore, x^^'^^ is on the tangent X2 = —2exi, below x~ . Similarly, x^^''^^ is below 
x+. We conclude that x^^''^^ G f^e. 

We have completed the proofs of the following sequence of lemmas. 

Lemma 7.9. The function (7.27), with fi± given by (7.26) for u = 0, is an optimizer for 
F+(x;-oo,0) — > To{x) i — F-(x;0,oo). 
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Lemma 7.10. The function (7.28), with u = e^^^l^ and fi± given by (7.26) for u = ^, is an 

optimizer for Lq{x) — > e, ^) — > T{x; ^) < — F~{x; ^, oo). 



Lemma 7.11. The function (7.29), with fi± given by (7.26) for u = 0, is an optimizer for 
L+{x;-£) — > To{x) i — L'{x;e). 

Putting together Lemmas |7.3[]7. 11 we obtain 

Theorem 7.12. For any x G 0^, we have 

Forp>2: Ne,p{x) < B,^p{x), Me,p{x) > 6,,p(x), 

For 1 < p < 2 : M,Jx) < B,^p{x), N,^p{x) > 
For 0<p<l: Pe,p{x) < B,^p{x), Re,p{x) > 

8. Proofs of the main inequalities 
Theorems 6.5 and 7.12 give us the exphcit expressions for Be,p{x) and he^p{x) for all p > : 
Theorem 8.1. For any x G 17^, we have 

Forp>2: Be,p{x) = N.^pix), b^^pix) = Ms,p{x), 

Forl<p<2: B,^p{x) = Me,p{x), be,p{x) = N^^pix), 
For <p <1 : Bs,p{x) = Pe,p{x), be,p(x) = Re,p{x). 

We are now in a position to prove all the theorems stated in Section [2] 

Proof of Theorem\2.1\ It suffices to consider Q = (0,1). Take any ip G BMO(Q) and let 



e = IIv'IIbmo(Q), xi = X2 = )q- Then 

be,p{xi,X2) < {\^\^)q < B^^p{xi,X2). 

Replacing ip with p — {p)^) gives 



^.11 



&,,p(0, X2-xl)<{\p- {p)^ \P)^ < S,,p(0, X2 



We now invoke Theorem |8.1[ for which we need the exact expressions for candidates M, N, P, 
and R. They come, respectively, from (5.15), (5.18), (5.27), and ( 5.30| ): 



Me,p{0,X2) = X^/'; iVe,p(0,X2) = ^T{p)eP-^X2; Pe,p(0,X2) 



; Re,piO,X2) = 2P-hP-^X2. 



Plugging these into (8.1) yields the stated inequalities. Furthermore, these inequalities are 
sharp because each becomes an equality, if we take ip to be the corresponding optimizer (/7(o,e2) 
from Section [7j Specifically, let 



"Pi 



'-e; te (o,i), 

e; te{l,l), 



elog(4t), tG{0,l), 
^-elog(4-4t), te (f,l). 



'-2e, te (0,|), 
¥'2 = <;u, ^3=lo, 

2e, tG(i,l). 

These are the optimizers given by (7.3), ( 7.27 ), and (7.29), respectively, each constructed for 
the point (0, e^). As was shown in Section 7, ||<y5j!c||BMO(Q) = (^'^^)q ~ ^'^^q-'^''^ ~ ^' ^ = 1> 2, 3. 
On the other hand, 

(l^i-(^i)ol% = (l^il%=^^ 



and 



(|V^2-(^2)q|% = (|^2|%=2.^. 

(1^3 -(9^3>J% = (19^31% = 



logt\Pdt= -T{p+l)eP, 



2 • - . {2ef = 2P-^eP. 



□ 
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Proof of Theorem 2^. Again, set Q = (0, 1). Take any (p G BMO(Q) and let e = ||9?||bmo(Q)- 

}j, xi = (^2)^. 



For any subinterval J oi Q let xf = (cp) j, X2 = {'p'^) j- Arguing as in the previous proof, we 



have 

be,p{0,xi - {xif) < {\^ - j < BeA^,xi - [xif). 

Let us consider these inequalities separately. For the one on the right, using that Se^p(0, 
is increasing, we have 

{\^-{^)jy)j<B,^,{Q,e^). 
Taking the supremum over all J, we get 



For each p this inequality is sharp: for p < 2 it is attained for Lp = ipi from (8.2), and for 
p > 2 it is attained for ip = p2- 

For the inequality on the left, take a sequence { Jn} of subintervals of Q such that 



lim (X2" - (a^i )^) 



n— >oo 



Using the continuity of be.p(0, •), we have 

6e,p(0,e2) = Jkn be,p(0,4" - {xi^f) < li^supd^^ - {^) jj) < II</'IIbmo''(Q)- 



This inequality is sharp for p > 2 : it is attained, again, for ip = ipi from (8.2). □ 

Remark 8.2. The last calculation in the proof shows why (p2 and p}^ cannot be used to show 
sharpness of the norm estimates in the cases 1 < p < 2 and < p < 1, respectively. Consider 
p = 1. The issue is that, while each function attains its BMO^ norm on (0, 1), neither attains 
its BMO^ norm on this interval. Indeed, one can easily calculate that both functions have 
1 -oscillations equal to e/2 on (0,1). However, 

Proof of Theorem \2.4\ We can set Q = (0,1). If ||9?||bmo(Q) = 0; there is nothing to prove. 
Assuming this is not the case, we use Theorem |2.1| to get 



2llwll^"P2 2llo3ll^"^i 



BMO(Q),| , > |p2> ^ / 2\ / \2 ^ ll^llBMO(Q),| , > .p^s 



Each of these inequalities becomes an equality when ip = ip2 from ( |8.2[ ). The left-hand side 
inequality in the statement of the theorem is attained, for instance, for ip = pi. □ 



Proof of Theorem \2.5[ Set Q = (0,1); all averages will be over Q. Let e = \\p\\bmo(Q) ^-iid 
assume e < 1. We have 

00 ^ 00 ^ 

^.Iv^-MI) = -dy, - {p)\') < 1 + B,,i(0,e2) + ^B,,,(0,e2) 

A:=0 ■ k=2 

00 00 1 — ^ 

= 1 + M,,i(0, s') + ^ yN,,,{0, e^) = l + e + -Y,e'' = 

k=2 ' k=2 



On the other hand, taking p = p2 from (8.2), we get 



/glv2-(^2>|^ 2-- t e'\'°^'\dt + 2- ^ - ^ ^ 



4 l-e 

This calculation also shows that the bound eo = 1 is sharp. □ 
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9. Other choices of / 

Throughout the paper, we have concentrated on one specific boundary function / : f{s) = 
\s\P, p > 0. However, the machinery developed in these pages works for many other choices of 
/. Let us briefly describe several such choices and their Bellman functions without going into 
details. 



9.1. /(s)=log|s| 



As mentioned earlier, this function corresponds to the case p 
lim{\ip\P)^/P = e<i°sl^l>. 

p-s>0 



0, since 



log|s|. Thus, F"*" and 



To 



It is easy to show that the corresponding Bellman functions are 
(9.1) Bsfl{x) = Pe,o{x) and bsfl{x) = -oo. 

Here P^^ is given by (5.26), with every block re-specified for /(s) 

F~ are given by (5.1) and (5.7), respectively; T is given by (4.27); and Lq, given by_(4.13) 
is simply ^logX2. To show that this is a viable global candidate, one needs Lemma 
prove the statement for B, use the local concavity of the candidate Psfi to run the induction 
of Section |6] and then apply the optimizer for P^^p from Section [7j To prove the statement for 
b, simply use the optimizer for R^ p from Section [?[ 
From ( |9.1[ ), we have sharp inequalities for ip G BMOg : 

- oo < {log\^\) < Pe,oi{^),{ipf), 



(9.2) 
and so 
(9.3) 



< e(iog|^-(^>|> < gP.,o(o,(^2^-M2) ^ gilog((v2^-M2) < 



While the second inequality in (9.2) is non-trivial and gives a sharp estimate on (log|9?|) for 
any pair of specified averages of if, its immediate consequence, the second-from-left inequality 



in (9.3), is not interesting, as it simply expresses the norm monotonicity (1.1). However, the 
leftmost inequality in (9.3) is important, as its sharpness means that BMO C BMO°. It is also 
the same result as one gets from the top line in Theorem 2.2, by taking the limit as p — )• 0"''. 

\s\P, p < 0. For this case, the situation reverses, compared to the previous one. 



9.2. f{s) 
We have 



Be,p{x) = OO, 



with P^^p given by (5.26) and so (5.27), and so 

PeMv)Av?)<M'') 



< CO, 



which produces the sharp inequalities 



9.3. fis) 



\s\, e < 1. This is the function that implicitly allowed us to prove the 



John-Nirenberg estimates of Theorem |2.5[ As we saw in Section [8} the key fact is that all 
upper Bellman functions B^^p and, separately, all lower Bellman functions b^^p have identical 
optimizers for p > 2. Since the Taylor expansion for e'*' — |s| is missing the term corresponding 
to p = 1, we expect the Bellman foliation of Qe to be the same as for \s\p, p > 2. Indeed, we 
can easily show that 

(9.4) B,j{x) = N,Ax), b,j{x) = M,j{x), 



where Mgj and A'gj are given by (5.14) and (5.17), with their blocks re-specified for this 
choice of /. Therefore, after a small bit of calculation, we have the sharp inequalities 



,vm^ - A¥^)^W < (ei^-<^>i> - (1^ - (^>i> < ^ 



2(1 



+ 1 < - 



e + 
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In fact, if / is any linear combination of powers greater than or equal to 2, with non- negative 
coefficients, the Bellman functions will be given by (9.4). If, on the other hand, / is a linear 
combination of powers between 1 and 2 with non-negative coefficients, the upper and lower 
Bellman functions will switch, i.e. they will be given by M^j and N^j, respectively. Further 
generalizations along these lines are possible. 
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